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FOREWORD 


This  report  presents  results  of  work  conducted  by  the  Boeing  Company,  Seattle, 
Washington,  under  Air  Force  Contract  F33615-77-C-3054  "Application  of  the  EASY 
Dynamic  Program  to  the  Analysis  of  Air  Cushion  Systems  on  Aircraft"  during  the 
period  from  15  April  1977  to  1  June  1979.  This  contract  was  conducted  under 
the  sponsorship  of  the  Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson 
Air  Force  Base,  Ohio  with  Peter  Skele  and  Lt.  David  L.  Fischer  as  project 
engineers. 

This  report  is  comprised  of  three  volumes. 

Volume  I  -  Component  Mathematical  Models 

Volume  II  -  Component  Computer  Programs 

(Parts  I  &  II) 

Volume  III  -  Description  of  Simulations 

In  addition,  a  User's  Manual  (Reference  1)  has  been  written  to  provide  a 
concise  reference  for  day  to  day  usage. 

The  results  presented  were  developed  by  the  Boeing  Aerospace  Company.  The 
program  managers  were  A.  J.  P.  Lloyd,  H.  H.  Straub  and  J.  R.  Kilner.  The 
principal  investigators  were  M.  K.  Wahi ,  G.  S.  Duleba,  J.  R.  Kilner  and 
P.  R.  Perkins. 
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coefficient  with  respect  to  r 

The  yawing  moment  coefficient  produced  by  the  ACLS  hardware 

The  nondimensional  derivative  of  the  aerodynamic  yawing  moment 
coefficient  with  respect  to  the  sideslip  angle. 

The  nondimensional  derivative  of  the  aerodynamic  yawing  moment 
coefficient  with  respect  to  B 

The  aerodynamic  yawing  moment  rudder  effectiveness  factor 

The  nondimensional  derivative  of  the  aerodynamic  yawing  moment 
coefficient  with  respect  to  the  aileron  deflection 
The  nondimensional  derivative  of  the  aerodynamic  yawing  moment 
coefficient  with  respect  to  the  rudder  deflection 
specific  heat  at  constant  pressure 

compressor  pressure  ratio 

The  coefficient  of  drag  produced  by  the  ground  effects 

The  nondimensional  X  axis  component  of  the  aerodynamic  bias 
force  coefficient 

The  nondimensional  derivative  of  the  X  axis  component  of  the 

aerodynamic  force  coefficient  with  respect  to  the  spoiler 
defl ectlon 

The  coefficient  of  drag  produced  by  the  ACLS  hardware 
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The  nondimensional  derivative  of  the  X  axis  component  of  the 
aerodynamic  force  coefficient  with  respect  to  u 
The  nondimensional  derivative  of  the  X  axis  component  of  the 
aerodynamic  force  coefficient  with  respect  to  the  angle  of  attack 
The  nondimensional  derivative  of  the  X  axis  component  of  the 
aerodynamic  force  coefficient  with  respect  to  the  elevator 
deflection 

The  nondimensional  derivative  of  the  Y  axis  component  of  the 

aerodynamic  force  coefficient  with  respect  to  p 

The  nondimensional  derivative  of  the  Y  axis  component  of  the 

aerodynamic  force  coefficient  with  respect  to  r 

The  nondimensional  derivative  of  the  Y  axis  component  of  the 

aerodynamic  force  coefficient  with  respect  to  the  spoiler 

deflection 

The  coefficient  of  sideforce  produced  by  the  ACLS  hardware 

The  nondimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  coefficient  with  respect  to  the  sideslip  angle 
The  nondimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  coefficient  with  respect  to  8 
The  sideforce  rudder  effectiveness  factor 

The  nondimensional  derivative  of  the  Y  axis  component  of  the 

aerodynamic  force  coefficient  with  respect  to  the  aileron 
defl ection 

The  nondimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  coefficient  with  respect  to  the  rudder 
defl ection 

The  coefficient  of  lift  produced  by  the  ground  effects 

Th«  nondimensional  Z  axis  component  of  the  aerodynamic  bias 
force  coefficient 
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DL 


The  nondimensional  derivative  of  the  Z  axis  component  of  the 
aerodynamic  force  coefficient  with  respect  to  q 
The  nondimensional  derivative  of  the  Z  axis  component  of  the 
aerodynamic  force  coefficient  with  respect  to  the  spoiler 
defl ection 

The  coefficient  of  lift  produced  by  the  ACLS  hardware 

The  nondimensional  derivative  of  the  Z  axis  component  of  the 

aerodynamic  force  coefficient  with  respect  to  u 

The  nondimensional  derivative  of  the  Z  axis  component  of  the 

aerodynamic  force  coefficient  with  respect  to  the  angle  of  attack 

The  nondimensional  derivative  of  the  Z  axis  component  of  the 

aerodynamic  force  coefficient  with  respect  to  a 

The  nondimensional  derivative  of  the  Z  axis  component  of  the 

aerodynamic  force  coefficient  with  respect  to  the  elevator 
defl ection 

The  nondimensional  derivative  of  the  Z  axis  component  of  the 
aerodynamic  force  coefficient  with  respect  to  the  stabilizer 
defl ection 

Horizontal  distance  between  center  of  free  shape  r ^  and  inboard 
attach  point 

Vertical  distance  between  center  of  free  shape  r ^  and  inboard 

attach  point 

Array  of  element  width 

Length  of  bump,  feet,  for  AM00E=1 

Incremental  distance  between  points  for  AM0DE=2 

Hydraulic  energy  absorber  damping  coefficient 

Air  bag  damping  coefficient  as  a  function  of  air  bag  flattened 

area 

Trunk  damping  coefficient  as  a  function  of  trunk  flattened  area 

Engine  compressor  temperature  rise 

Area  adjustment  factor  for  perforation  area  AHA1 
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Area  adjustment  factor  for  perforation  area  AHA2 
Cable  modulus  of  elasticity 

Area  adjustment  factor  for  perforation  area  AHC1 
Area  adjustment  factor  for  perforation  area  AHC2 
Corrected  engine  Speed 
Actual  engine  pressure  ratio  (Forward) 

Connand  engine  pressure  ratio 

Actual  engine  pressure  ratio 

Actual  engine  pressure  ratio  (reverse) 

Tape  stress-strain  function 
Element  cushion  pressure  force 
Element  damping  force 
Element  friction  force  in  X  axis 
Element  friction  force  in  Y  axis 

Table  look  up  function  for  linear  interpolation  on  one 
independent  variable 

Table  look  up  function  for  curvilinear  to  rectilinear 
transformation  and  linear  interpolation  on  two  independent 
variables 

Cable  impact  load,  left  side 
Cable  impact  load,  right  side 
Element  ground  contact  reaction 
X  axis  force  summation  term 

The  X  stability  axes  or  body  axes  component  of  the  resultant 
aerodynamic  force 

X  body  axes  component  of  the  sum  of  all  the  external  forces 
acting  on  the  aircraft 
Y  axis  force  summation  term 

The  Y  stability  axes  or  body  axes  component  of  the  resultant 
aerodynamic  force 
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FY1 

Fz 

FZaero 

FZ1 

9 


iDR 

IM 


IT 


xy 

*xz 

'y 

‘yz 

'z 

k 


'GE 


L 

KCl 

KCM, 


Y  body  axes  component  of  the  sum  of  all  the  external  forces 
acting  on  the  aircraft 
Z  axis  force  summation  term 

The  Z  stability  axes  or  body  axes  component  of  the  resultant 
aerodynamic  force 

Z  body  axes  component  of  the  sum  of  all  the  external  forces 

acting  on  the  aircraft 

Acceleration  of  gravity 

Initial  cable  height  above  the  runway 

Hook  pivot  to  vehicle  c.g.  vertical  distance 

Hook  to  vehicle  c.g.  vertical  distance 

Inertia  of  tape  stored  on  drum 

Drum  inertia 


Logical  Switch 


Logical  Switch 


1  =  Hook  mised  cable 
0  =  Hook  engaged  cable 

1  =  Down  range  of  arresting  system 
0  =  Up  range  of  arresting  system 


X  body  axes  moment  of  inertia 

X-Y  body  axes  cross  products  of  inertia 

X-Z  body  axes  cross  products  of  inertia 

Y  body  axes  moment  of  inertia 

Y-Z  body  axes  cross  products  of  inertia 

Z  body  axes  moment  of  inertia 

Polytropic  constant  for  air 

The  ground  effects  factor 


Rolling  moment  factor 


The  rolling  moments  sideslip  factor 
The  pitching  moment  sideslip  factor 


xlx 


LIST  OF  SYMBOLS,  SUBSCRIPTS  AND  SUPERSCRIPTS  (CONTINUED) 


SYMBOLS 


*8 


KCX 
K  r 


KCV 


t8 


KCi 


•8 


KD 

K, 

lb 


'M 


8 


’8 


*8 

re 


:e 


Yawing  moment  factor 

The  yawing  moment  sideslip  factor 

The  X  axis  sideslip  factor 

Si  deforce  factor 

The  Y  axis  sideslip  factor 

The  Z  axis  sideslip  factor 

Spin  down  gain  over  the  spin  down  time  constant 
The  ground  effects  factor 

The  rolling  moment  sideslip  factor 
The  pitching  moment  sideslip  factor 
The  yawing  moment  sideslip  factor 
The  X  axis  sideslip  factor 
The  Y  axis  sideslip  factor 
The  Z  axis  sideslip  factor 

Aerodynamic  rolling  moment 
Characteristic  length 
Membrane  meridian  length 
Membrane  meridian  free  length 
Rolling  moment  factor  for  ground  effects 
Hook  arm  length 

Unflattened  length  segment  on  outboard  side  contributing  flow  to 
atmosphere 
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L 


ha  I 


-hb 


Unflattened  length  segment  contributing  flow  to  atmosphere 
Unflattened  length  segment  on  inboard  side  contributing  flow  to 
atmosphere 

Unflattened  length  segment  contributing  flow  to  cushion 
Element  width  of  perforated  area  array 
Element  membrane  peripheral  length 
Membrane  merdian  installed  length 

The  dimensional  derivative  of  the  aerodynamic  rolling  moment 
with  respect  to  p 

Element  peripheral  distance  from  outboard  attach  point  (B)  to 
edge  of  perforated  area 

The  dimensional  derivative  of  the  aerodynamic  rolling  moment 
with  respect  to  r 

The  dimensional  derivative  of  the  aerodynamic  rolling  moment 


with  respect  to  the  spoiler  deflection 

The  rolling  moment  coefficient  produced  by  the  ACLS  hardware 


The  X  body  axes 
The  Y  body  axes 
The  Z  body  axes 
Length  of  trunk 
Length  of  trunk 
Length  of  trunk 
Length  of  trunk 
Length  of  trunk 
Length  of  trunk 
The  dimensional 
with  respect  to 
The  dimensional 
with  respect  to 


gust  wind  scale 
gust  wind  scale 
gust  wind  scale 

or  air  bag  segment  inscribed  by  angle*  ^ 
or  air  bag  segment  inscribed  by  angle*  ^ 
material  crumpled  against  the  ground 
or  air  bag  element  footprint 

element  footprint  contributing  flow  to  atmosphere 
element  footprint  contributing  flow  to  cushion 
derivative  of  the  aerodynamic  rolling  moment 
the  sideslip  angle 

derivative  of  the  aerodynamic  rolling  moment 

0 


The  dimensional  derivative  of  the  aerodynamic  rolling  moment 
with  respect  to  the  aileron  deflection 
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M 

M 

m 

MF 

M, 


GE 


M_ 


’SP 


TR 


’B 


FS 


'GE 


The  dimensional  derivative  of  the  aerodynamic  rolling  moment 
with  respect  to  the  rudder  deflection 

Aerodynamic  pitching  moment 

Mach  number 

Aircraft  mass 

Dimensionless  shear  wind  factor 


The  pitching  moment  produced  by  the  ground  effects 
The  bias  pitching  moment 

The  dimensional  derivative  of  the  aerodynamic  pitching  moment 
with  respect  to  q 

The  dimensional  derivative  of  the  aerodynamic  pitching  moment 

with  respect  to  the  spoiler  deflection 

The  pitching  moment  produced  by  the  ACLS  hardware 

The  dimensional  derivative  of  the  aerodynamic  pitching  moment 

with  respect  to  u 

The  dimensional  derivative  of  the  aerodynamic  pitching  moment 
with  respect  to  the  angle  of  attack 

The  dimensional  derivative  of  the  aerodynamic  pitching  moment 
with  respect  to  a 

The  aerodynamic  pitching  moment  for  sideslip 


The  dimensional  derivative  of  the  aerodynamic  pitching  moment 
with  respect  to  the  elevator  deflection 

The  dimensional  derivative  of  tne  aerodynamic  pitching  moment 
with  respect  to  the  stabilizer  deflection 
Aerodynamic  yawing  moment 

The  dimensional  derivative  of  the  aerodynamic  yawing  moment  with 
respect  to  the  spoiler  deflection 
Yawing  moment  factor  for  ground  effects 

The  dimensional  derivative  of  the  aerodynamic  yawing  moment  with 
respect  to  p 
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'TR 


B 


N 


«♦ 

N. 


amb 

Pc 

PPU 

pr 

Pram 

PRAT 

PT 

Pt 

Pt1 

Ptr 

P2 


The  dimensional  derivative  of  the 
respect  to  r 

The  yawing  moment  produced  by  the 
The  dimensional  derivative  of  the 
respect  to  the  sideslip  angle 
The  dimensional  derivative  of  the 
respect  to  6 

The  dimensional  derivative  of  the 
respect  to  the  aileron  deflection 
The  dimensional  derivative  of  the 
respect  to  the  rudder  deflection 
Hoop  load 

Combined  hoop  load 

Merdian  load 
Combined  meridian  load 


aerodynamic  yawing  moment  with 
ACLS  hardware 

aerodynamic  yawing  moment  with 
aerodynamic  yawing  moment  with 

aerodynamic  yawing  moment  with 
aerodynamic  yawing  moment  with 


X  body  axes  component  of  the  angular  velocity  of  the  aircraft 

c.g.  with  respect  to  the  earth  axes 

Ambient  air  pressure 

Ambient  air  pressure 

Cushion  air  pressure 

Engine  bleed  port  inlet  pressure 

Cushion  to  trunk  pressure  ratio  (pc-f>a)/(Pt-Pa) 

Ram  pressure 

Engine  port  pressure  ratio 
Ram  pressure  at  engine  inlet 
Trunk  or  air  bag  air  pressure 
Left  bag  air  pressure 
Right  bag  air  pressure 
Engine  exhaust  pressure 
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SYMBOLS 


4a2 

QC2 


R 

r 


S 

s 

Si 

SPD 

T 

t 

^amb 

TCO 

Tcu 

TF 

TH 

THR 

T. 


ram 


Y  body  axes  component  of  the  angular  velocity  of  the  aircraft 

c.g.  with  respect  to  the  earth  axes 

Total  flow  from  the  flattened  area  to  the  atmosphere 

Total  flow  from  the  flattened  area  to  the  cushion 

Gas  constant 

Z  body  axes  componeni.  of  the  angular  velocity  of  the  aircraft 
c.g.  with  respect  to  the  earth  axes 
Perforation  flow  area  adjustment  parameter 
Drum  hub  radius 

Instantaneous  radius  of  stored  tape,  left  side 
Perforation  flow  area  adjustment  parameter 
Drum  outside  radius 

Instantaneous  radius  of  stored  tape,  right  side 

Radius  of  curvature  for  segment  Lj 

Radius  of  curvature  for  segnent  L2 

Wing  area 

Laplace  operator 

User  supplied  element  scaling  factors  array 
Engine  speed 

Time  of  events  (simulation) 

Tape  thickness 

Ambient  temperature 

Engine  time  constant 

Temperature  of  air  flow  into  cushion 

Forward  engine  thrust 

thrust 

Required  thrust  level 
Tape/cable  tensile  force,  left  side 
Tape/cable  tensile  force,  right  side 
Reverse  engine  thrust 
Ram  temperature 
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SYMBOLS 


TT 


'  tl 

Ttr 

Ttr 

Tx 

TX1 


TX 


aero 


Ty 

TYl 


TY 


aero 


Tz 

^aero 

TZ1 


T2 

u 


UUS 

V 

7 

v 


Ram  temperature  at  engine  inlet 

Tension  in  trunk  material,  lbs  in  tangential  direction  per  foot 

length  in  the  axial  direction 

Temperature  of  air  flow  into  left  air  bag 

Temperature  of  air  flow  into  trunk 

Temperature  of  air  flow  into  right  air  bag 

Roll  moment  sumnation  term 

X  body  axes  component  of  the  sum  of  all  the  external  torques 
acting  about  the  aircraft  c.g. 

The  X  stability  axes  or  body  axes  component  of  the  resultant 

aerodynamic  torque 

Pitch  moment  sumnation  term 

Y  body  axes  component  of  the  sum  of  all  the  external  torques 
acting  about  the  aircraft  c.g. 

The  Y  stability  axes  or  body  axes  component  of  the  resultant 

aerodynamic  torque 

Yaw  moment  sumnation  term 

The  Z  stability  axes  or  body  axes  component  of  the  resultant 
aerodynamic  torque 

Z  body  axes  component  of  the  sum  of  all  the  external  torques 
acting  about  the  aircraft  c.g. 

Engine  exhaust  temperature 

X  body  axes  component  of  the  linear  velocity  of  the  aircraft 
c.g.  with  respect  to  either  the  earth  axes  or  the  atmosphere 
The  X  body  axes  component  of  the  wind  shear  vector 
Airspeed 

The  velocity  of  the  aircraft  c.g.  with  respect  to  the  atmosphere 

Y  body  axes  component  of  the  linear  velocity  of  the  aircraft 
c.g.  with  respect  to  either  the  earth  axes  or  the  atmosohere 
Cushion  volume 
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'tl 


'tr 

vtl* 


*tr* 

VWS 

w 

w 


al 

war 
WB1  eed 
Wca 

Wcu 

WLCG 

WLT 

wn 

Wta 

Mtc 

Wt1 

Wtr 

wws 

W2 

X 


Xa1 

xBai 


Calculated  cushion  volume 

Hook  velocity  in  Earth  axes 
Trunk  volume 
Left  air  bag  volume 
Right  air  bag  volume 
Calculated  left  air  bag  volume 
Calculated  right  air  bag  volume 

Calculated  trunk  volume 

The  Y  body  axes  component  of  the  wind  shear  vector 
Tape  width 

Z  body  axes  component  of  the  linear  velocity  of  the  aircraft 

c.g.  with  respect  to  either  the  earth  axes  or  the  atmosphere 

Air  flow  rate,  left  air  bag  to  atmosphere 

Air  flow  rate,  right  air  bag  to  atmosphere 

Engine  bleed  air  flow  rate 

Air  flow  rate  cushion  to  atmosphere 

Air  flow  rate  supply  to  cushion 

Water  line  center  of  gravity  location 

Water  line  trunk  axes  location 

natural  frequency 

Air  flow  rate  trunk  to  atmosphere 

Air  flow  rate  trunk  to  cushion 

Air  flow  rate  supply  left  air  bag 

Air  flow  rate  suppply  to  trunk,  or  right  air  bag 

The  Z  body  axes  component  of  the  wind  shear  vector 

Engine  exhaust  air  flow  rate 

The  X  earth  axes  component  of  the  vector  from  the  origin  of  the 
earth  axes  to  the  aircraft  c.g.  (i.e.  position  vector) 

Array  of  X  coordinates  of  the  element  inboard  attachment  points 
An  X  axis  term  relating  XA  and  XBT 
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A  term  relating  YO,  RA  and  BET 
X  body  axes  position  of  point  "t" 

X  earth  axes  position  of  point  "t" 

The  X  body  axes  component  of  the  velocity  of  point  "t"  relative 
to  the  earth  axes 

Drag  produced  by  the  ground  effects 

The  dimensional  X  component  of  the  aerodynamic  bias  force 
The  dimensional  derivative  of  the  X  axis  component  of  the 
aerodynamic  force  with  respect  to  the  spoiler  deflection 
Drag  produced  by  the  ACLS  hardware 

The  dimensional  derivative  of  the  X  axis  component  of  the 
aerodynamic  force  with  respect  to  u 

The  dimensional  derivative  of  the  X  axis  component  of  the 
aerodynamic  force  with  respect  to  the  angle  of  attack 

The  dimensional  derivative  of  the  X  axis  component  of  the 
aerodynamic  force  with  respect  to  the  elevator  deflection 
The  Y  earth  axes  component  of  the  vector  from  the  origin  of  the 
earth  axes  to  the  aircraft  c.g.  (i.e.  position  vector) 

Array  of  Y  coordinates  of  the  element  inboard  attachment  points 
An  Y  axis  term  relating  YA  and  YBT 
A  term  relating  YBA  and  YO 
Y  body  axes  position  of  point  "t" 


The  Y  body  axes  component  of  the  velocity  of  point  "t"  relative 
to  the  earth  axes 

SI  deforce  factor  for  ground  effects 
Spanwise  intercept,  left  side 
Tape  drum  to  sheave  distance 

Horizontal  distance  between  trunk  or  air  bag  center  line  and 
Inboard  attach  point 
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^BT 

zEg 

ZEt 


Free  shape  Yo 

The  dimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  with  respect  to  p 
Spanwise  intercept,  right  side 

The  dimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  with  respect  to  r 
Runway  span  between  sheaves 

The  dimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  with  respect  to  the  spoiler  deflection 
The  si  deforce  produced  by  the  ACLS  hardware 
The  dimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  with  respect  to  the  sideslip  angle 
The  dimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  with  respect  to  0 

The  dimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  with  respect  to  the  aileron  deflection 

The  dimensional  derivative  of  the  Y  axis  component  of  the 
aerodynamic  force  with  respect  to  the  rudder  deflection 

The  Z  earth  axes  component  of  the  vector  from  the  origin  of  the 
earth  axes  to  the  aircraft  c.g.  (i.e.  position  vector) 

Array  of  Z  coordinates  of  the  element  inboard  attachment  points 
A  Z  axis  term  relating  ZA  and  Z8T 

Arbitrary  constant  in  body  axes 
Z  body  axes  position  of  point  "t" 

Terrain  elevation  at  point  "t" 

Z  earth  axes  position  of  point  "t" 

The  Z  body  axes  component  of  the  velocity  of  point  "t"  relative 
to  the  earth  axes 

Trunk  or  air  bag  ground  gap  height 


xxvl 1 1 


LIST  OF  SYMBOLS,  SUBSCRIPTS  AND  SUPERSCRIPTS  (CONTINUED) 


SYMBOLS 


za 


The  lift  produced  by  ground  effects 

The  dimensional  Z  axis  component  of  the  aerodynamic  bias  force 
Vertical  distance  between  the  aircraft  hard  structure  and  bottom 
of  the  trunk  or  air  bag 

Free  shape  ZO,  the  point  at  which  the  membrane  first  makes 
contact  with  the  ground  for  a  given  Pr 
Unsealed  value  of  ZOFS 
Unsealed  value  of  ZO 

The  dimensional  derivative  of  the  Z  axis  component  of  the 
aerodynamic  force  with  respect  to  q 

The  dimensional  derivative  of  the  Z  axis  component  o*  the 
aerodynamic  force  with  respect  to  the  spoiler  deflection 
The  lift  produced  by  the  ACLS  hardware 
The  dimensional  derivative  of  the  Z  axis  component  of  the 
aerodynamic  force  with  respect  to  u 

The  dimensional  derivative  of  the  Z  axis  component  of  the 
aerodynamic  force  with  respect  to  the  angle  of  attack 
The  dimensional  derivative  of  the  Z  axis  component  of  the 
aerodynamic  force  with  respect  to  d 
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SECTION  I 


INTRODUCTION 


The  Air  Force  Flight  Dynamics  Laboratory  (AFFOL)  has  been  involved,  for  a 
number  of  years,  in  developing  Air  Cushion  Landing  Systems  (ACLS).  The  ACLS, 
though  not  completely  developed  into  a  mature  technology,  holds  great  promise 
as  a  viable  alternative  to  the  traditional  wheel  landing  gear  system. 

An  essential  tool  for  the  development  of  a  successful  ACLS  is  a  general 
purpose  computer  program  for  dynamic  simulation  and  analysis.  To  provide  this 
dynamic  simulation  and  analysis  capability  the  Boeing  Military  Airplane 
Development  (BMAD)  Division  of  the  Boeing  Aerospace  Company  was  awarded  a 
contract  to  create  an  ACLS  component  library  for  the  EASY  Model  Generation  and 
Analysis  Computer  program. 

The  EASY  computer  program  is  an  analytical  tool  which  has  been  constructed  to 
assist  in  the  development  of  improved  control  systems.  The  program  has  many 
capabilities  which,  until  now,  have  not  been  available  for  dynamic  analysis  or 
have  been  available  only  in  separate  programs.  Capabilities  of  the  program 
Include: 

o  Steady  state  analysis 
o  Non-linear  simulation  (time  histories) 
o  Linear  model  generation 
o  Stability  matrix  calculation 

o  Eigenvalue  and  eigenvalue  sensitivity  determination 
o  Stability  margins 

o  Transfer  functions  (Bode,  Nyquist,  and  Nichols  plots) 
o  Root  locus 
o  Algebraic  sensitivity 
o  Optimal  control  synthesis 
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The  major  tasks  of  the  EASY-ACLS  contract  are  as  follows: 


(1)  Develop  Aerodynamic  Library  Components 

(2)  Simulate  Jindivik  Flight 

(3)  Develop  ACLS  Library  Components 

(4)  Simulate  Jindivik  Drop  Test 

(5)  Simulate  Jindivik  Takeoff  and  Landing 

(6)  Simulate  XC-3A  Flight 

(7)  Simulate  XC-3A  Landing,  Takeoff  and  Taxi 

(3)  Simulate  YC-14  Flight 

Task  1  of  the  contract  required  the  development  of  aircraft  mathematical 
modeling  components  consisting  of  Equations  of  Motion,  Wind  Modeling, 
Aerodynamic  Force  and  Moment  Modeling,  and  Engine  Modeling. 

Task  2  of  the  contract  required  modeling  the  Jindivik  RPV  in  free  flight 
(without  controls)  and  in  clean  configuration  (without  trunk)  using  the 
components  developed  in  Task  1;  and  to  determine  the  longitudinal  and  lateral 
stability  of  the  vehicle. 

Task  3  of  the  contract  required  development  of  EASY  components  required  for 
modeling  air  cushion  systems,  including  air  flow  system  components,  fan  and 
ejector  components,  inelastic  and  elastic  trunk  and  cushion  components, 
terrain  model  component,  arresting  system  component,  and  air  bag  skid 
component . 

Task  4  of  the  contract  required  that  the  EASY  program  capabilities  be 
demonstrated  by  simulating  a  Jindivik  drop  test  using  the  models  developed  in 
Tasks  1  and  3.  It  further  required  that  the  simulation  results  be  correlated 
with  those  of  an  actual  drop  test  on  a  recovery  trunk  (ACRS  No.  2). 

Task  5  of  the  contract  required  that  the  EASY  program  capabilities  be  further 
demonstrated  by  landing  and  takeoff  simulations.  The  approach,  touchdown,  and 
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slideout  of  the  Jindivik  on  the  recovery  trunk;  along  with  the  roll  out, 
rotation  and  climbout  of  the  Jindivik  on  the  takeoff  trunk;  using  models 
developed  in  Tasks  1  and  3  were  simulated.  In  addition,  simulations  to  check 
out  the  airbag  skid  component  and  the  arresting  system  component  were  also 
performed. 

Task  6  of  the  contract  required  modeling  the  XC-8A  experimental  ACLS  aircraft 
during  trim  level  flight.  A  "saw  tooth"  pertubation  of  the  rudder  and 
elevator  were  made  to  evaluate  the  lateral  and  longitudinal  stability  of  the 
aircraft. 

Task  7  of  the  contract  required  that  the  EASY  program  capabilities  be 
demonstrated  by  simulating  the  XC-8A  aircraft  during  takeoff,  landing  and 
taxi . 

Task  8  of  the  contract  required  modeling  the  YC-14  aircraft  during  an  approach 
configuration.  A  square  impulse  deflection  of  the  rudder  and  elevator  were 
made  to  evaluate  the  lateral  and  longitudinal  stability  of  the  aircraft. 

The  resulting  EASY-ACLS  component  library  has  been  documented  in  three  main 
volumes.  Volume  I,  this  volume,  contains  the  mathematical  models  for  the 
standard  components  in  the  ACLS  library.  Tasks  1  and  3.  Volume  II  contains  a 
detailed  description  of  all  the  component  subroutines,  including  flowcharts 
and  detailed  listings.  Volume  III  contains  a  description  of  the  development, 
results  and  conclusions  of  EASY  computer  simulations  accomplished  during  this 
contract.  A  summary  document  (Reference  1)  has  also  been  prepared  as  a  User's 
Manual . 
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SECTION  II 


AIRCRAFT  MODELING  COMPONENTS 


Introduction 

The  mathematical  formulation  of  the  equations  of  motion  of  an  aircraft  through 
nonstationary  air  subject  to  aerodynamic,  gravitational  and  other  external 
forces  requires  two  sets  of  parametric  equations,  which  are: 

1.  The  force  equations  which  relate  the  linear  acceleration  of  the 
center  of  mass  to  the  external  forces. 

2.  The  moment  equations  which  relate  the  angular  acceleration  about  the 
center  of  mass  to  the  external  torques. 

These  two  sets  of  differential  equations  are  integrated  to  yield  both  the 
linear  and  angular  position  and  velocity  of  the  aircraft  center  of  mass. 

2.1  Axes  and  Euler  Angles 

The  coordinate  axes  systems  used  to  describe  the  motion  of  an  aircraft  are  the 
earth  fixed  axes,  body  axes,  vehicle  axes,  wind  axes  and  stability  axes. 

Since  these  are  standard  reference  frames  commonly  used  to  describe  aircraft 
motion  they  will  not  be  defined  here.  The  reader  who  is  unfamiliar  with  these 
coordinate  axes  systems  should  consult  Reference  2,  Chapter  4;  or  Reference  3. 

A  description  of  the  basic  aerodynamic  terms  used  in  the  equations  of  motion 
are  presented  in  the  succeeding  paragraphs. 

The  equations  of  motion  for  the  aircraft  are  expressed  as  parametric  equations 
resolved  into  the  body  axes  system.  Since  the  body  axes  (see  Figure  1)  are  a 
rotating  noninertial  reference  frame,  the  earth  axes  were  assumed  to  be  the 
inertial  reference  frame  for  the  derivation  of  the  equations  of  motion  of  the 
aircraft. 
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It  *  roiling  moment  p  =*  rate  of  roll 

if  —  pitching  moment  q  **  rate  of  pitch 

JST  »  yawing  moment  r  =*  rate  of  yaw 

[X,  T,  Z]  mm  component*  of  resultant  aerodynamic  force 
[a,  9,  w]  *  component*  of  velocity  of  C  relative  to  atmosphere 


SIGN  CONVENTION 

FOR 

POSITIVE  CONTROL 

SURFACE  DEFLECTIONS 

CONTROL  SURFACE 

+  DEFLECTION 

ELEVATOR 

DOWN 

STABILIZER 

LEADING  EDGE  UP 

SPOILER 

UP 

FLAPS 

DOWN 

RUDDER 

TO  THE  PILOTS  LEFT 

AILERON 

TIPS  RIGHT  WING  DOWN 

Figure  1  Notation  for  Body  Axes 


Figure  2  Plane  of  Synmetry  -  CXZ 
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To  derive  the  equations  of  motion  in  the  body  axes  requires  the  linear 
velocity  and  the  angular  velocity  of  the  body  axes  with  respect  to  the  earth 
axes.  The  symbols  which  are  commonly  used  to  represent  these  velocity  vectors 
are  given  below: 

(u,v,w)  *  X,  Y  and  Z  body  axes  components  of  the  linear  velocity  of  the 

aircraft  c.g.  with  respect  to  the  earth  axes. 

( p »q , r )  =>  X,  Y  and  Z  body  axes  components  of  the  angular  velocity  of 

the  c.g.  with  respect  to  the  earth  axes. 

Since  the  origin  of  the  body  axes  system  is  located  at  the  center  of  gravity 
(c.g.)  of  the  aircraft  these  terms  are  used  interchangably  throughout  the  text. 

The  position  of  the  aircraft  is  defined  as  the  vector  from  the  origin  of  the 
earth  axes  to  the  aircraft  c.g.  X,  Y,  and  Z  are  the  earth  axes  components  of 
this  position  vector.  The  aircraft  altitude  is  in  turn  defined  as: 

Altitude  =  -Z 

The  calculation  of  the  aerodynamic  forces  and  moments  present  in  the  equation 
of  motion  are  a  function  of  the  velocity  of  the  aircraft  c.g.  with  respect  to 
the  atmosphere  (V).  This  velocity  vector  is  commonly  resolved  into  body  axes 
components  (u,  v,  w).  The  airspeed  (V)  is  defined  as  the  magnitude  of  this 
velocity  vector. 

The  same  symbols  (u,  v,  w)  are  used  for  the  velocity  components  of  the 
aircraft  c.g.  with  respect  to  the  earth  axes.  These  velocity  components  are 
used  in  the  equations  of  motion  of  the  aircraft  (i.e.  Equation  2.2-1  and 
2.2-2)  whereas  the  components  of  V  are  used  in  evaluating  the  aerodynamic 
forces  and  moments  (see  Section  2.3).  It  should  be  noted  that  V  minus  the 
wind  vector  equals  the  velocity  of  the  aircraft  c.g.  with  respect  to  the  earth 
axes. 

The  orientation  of  the  stability  axes  and  the  wind  axes  systems  are  defined 
with  respect  to  the  airspeed  vector  V.  The  X  stability  axis  lies  along  the 
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projection  of  V  onto  the  plane  of  symmetry  of  the  aircraft  (See  Fig.  2).  The 
X  wind  axis  lies  along  the  V  velocity  vector.  The  relationship  of  the  wind 
axes  system  to  the  aerodynamic  lift  (L),  drag  (D)  and  sideforce  (C)  components 
of  the  total  aerodynamic  force  are  shown  in  Figures  3  and  4. 

The  aircraft  angle  of  attack  (a)  used  in  EASY  is  defined  as  the  angle  made  by 
the  X  stability  axis  and  the  X  body  axis.  The  sideslip  angle  (6)  is  defined 
as  the  angle  the  X  wind  axis  makes  with  the  plane  of  symmetry  of  the  aircraft 
(see  Figures  3  and  4)).  Since  (u,  v,  w)  are  the  X,  Y  and  2  body  axes 
components  of  V,  the  angle  of  attack  and  sideslip  are  given  by: 


a-  tan"1  ^  ,  where  -  a_<TT,  and 
B=  sin’1  j  ,  where  -  TK  g  £TT . 

If  V,  a  and  g  are  known,  then  u,  v  and  w  can  be  determined  from; 

u  -  VcosBcosa 
v  *  Vsing 
w  *  Vcosgsina 

The  flight  path  angle  (Y)  is  defined  as  the  angle  between  the  horizontal 
plane  and  the  X  wind  axis  (see  Figure  5). 


The  orientation  of  any  reference  frame  relative  to  another  can  be  given  by 
three  angles,  called  Euler  angles,  which  are  the  consecutive  rotations  about 
the  Z,  Y,  X  axes  that  carry  one  frame  into  coincidence  with  the  other.  In 
flight  dynamics,  the  Euler  angles  used  are  those  which  rotate  the 
vehicle-earned  vertical  frame  Fv  into  coincidence  with  either  the  body  axes 
or  the  wind  axes  systems.  The  Euler  angles  are  denoted  (4i,0,$)  for  body  axes, 
and  (il>w,  0W,  $w)  for  wind  axes.  Figure  6  shows  the  sequence  of  rotations. 


1.  A  rotation  about  Czv,  carrying  the  axes  to  ^  is  the 

azimuth  angle  (yaw). 


• t i. 


flifht  path 


Figure  6  The  Euler  Angles 


2 


.  A  rotation  0  about  Cy^,  carrying  the  axes  to  Cx^y^z^.  0  is  the 
elevation  angle  (pitch). 

3.  A  rotation  $  about  Cx^,  carrying  the  axes  to  their  final  position  Cxyz. 
♦  is  the  bank  angle  (roll). 

In  order  to  avoid  ambiguities  in  the  definition  of  the  angles  4*  ,9,  $  their 
ranges  are  limited  to: 

-TT<.4><  IT  Of  0  £ 4<<  2  IT 

<e<  | 

-7T£4>  <  IT  or  0  £<$  <  2  7T 

It  should  also  be  noted  that  the  Euler  angles  which  relate  the  wind  and  body 
axes  are  (  -0,  o ,  0) . 


2.2  Equations  of  Motion 


Assuming  the  earth  fixed  reference  frame  is  an  inertial  reference  frame  the 
following  six  parametric  non-linear  differential  equations  can  be  derived  to 
fully  describe  the  motion  of  a  vehicle.  For  a  detailed  derivation  of  these 
equations  the  reader  is  referred  to  Reference  2,  Chapter  4. 


Aircraft  Equations  of  Motion  Resolved  into  the  Body  Axis 


Force  Equations: 


FX1  +  Fxaero  cosa-FZaer()  sina-  mgsinO 
FY1  +  FYaer()  +  mgcos9sinjS 

FZ1  +  FZaero  cosa+pxaero  sina+  mgcosOcosd 


m(u+qw-rv) 
m(  v+ru-pw) 
m(w+pv-qu) 


2.2-1 
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Moment  Equations: 


TX1  +  TXaero  coso  -TZaero  sina*  lx  p  ~lxz  *  ^z)v 

TY1  +TYaero  -  Iy.q  -I„  (r2-P2)  -  (1,-I^rp  2.2-2 

TZ1  +  TZaero  cosa+TXaero  sin°a  lz  r  'lxz  {^r)  *  ^x'V 


where: 


FX1,  FY1  and  FZ1  are  the  sum  of  all  the  external  forces  resolved 
Into  the  body  axis, 

TX1,  TY1  and  TZ1  are  the  sum  of  all  the  extenal  torques  resolved  into  the 
body  axis. 


FX....,.,  FY, _  and  FZ.  are  the  components  of  resultant  aerodynamic 

aero’  aero  aero 

force  resolved  into  the  stability  axis, 

TX  _  TY  and  TZ  are  the  components  of  the  resultant  aerodynamic 

aero  aero  aero 

torque  resolved  into  the  stability  axis. 


m  -  aircraft  mass,  and 


g  ■  acceleration  of  gravity. 


The  equations  are  first  order  in  u,  v,  w,  p,  q,  and  r.  The  Euler  angles  and 
the  aircraft  position  vector  are  determined  by  the  following  kinematic 
relationships: 


i  •  p  +  q  sin$  tan  0  +  r  cos$  tan  0  (a) 

0  ■  q  cos$  -  r  sin$  (b) 

■  (q  sin$  +  r  cos$)sec0  (c) 

X  ■  ucosQcos'J' +  v(sin$sin0cos^  -  cos0sinij/)  (a) 

+  w(cos0sin0cos4»  +  sin0sin4O 
Y  ■  ucosOsin'J' +  v(sin0sin0simj>  +  cos^cos^)  (b) 

+  w(cos^sin0sini^  -  sin^cos^) 

Z  ■  -usin0+  vsin$cos0+  wcos$cos0  (c) 


2.2-3 


2.2-4 
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The  assumptions  used  in  the  derivation  of  these  equations  were: 

(1)  the  aircraft  is  a  rigid  body,  (2)  the  mass  of  the  aircraft  is  constant 
for  the  duration  of  the  analysis,  (3)  gravity  is  constant,  (4)  the  earth  is 
an  inertial  reference,  and  (5)  there  is  body  axis  synmetry  about  the  x-z  plane 
(i*e*.  1 xy  3  *yz  3  diagram  of  these  equations  is  given  in  Figure 

7. 


2.3  Aerodynamic  Forces  and  Moments 

In  general,  the  aerodynamic  forces  and  moments  present  in  the  equations  of 
motion  are  functions  of  air  density,  angle  of  attack,  the  time  derivatives  of 
the  angle  of  attack,  sideslip  angle,  the  time  derivative  of  the  sideslip 
angle,  linear  velocity,  the  time  derivative  of  the  linear  velocity,  angular 
velocity,  the  time  derivative  of  the  angular  velocity  and  the  control  surface 
deflections.  Since  the  aerodynamic  force  and  moment  equations  can  be  written 
in  either  dimensional  or  nondimensional  forms,  both  forms  are  presented  in  the 
succeeding  paragraphs. 

2.3.1  Dimensional  Aerodynamic  Force  and  Moment  Equations 

The  dimensional  aerodynamic  force  and  moment  equations  are  derived  by  assuming 
that  these  equations  are  continuous  functions  of  the  variables  listed  in  the 
preceeding  paragraph  and  that  the  aircraft  experiences  only  small 
perturbations  from  a  trirrmed  flight  configuration.  Based  on  these  assumptions 
a  Taylor  series  expansion  of  the  forces  and  moments  about  the  trimmed  flight 
configuration  can  be  written  as: 


Second  and  higher  order  terms  have  been  dropped  (i.e.  small  pertubations  were 
assumed).  The  terms  Fq  and  Mq  are  the  trimmed  flight  aerodynamic  forces  and 
moments.  The  first  order  derivative  terms  given  in  these  equations  are  the 
dimensional  aerodynamic  stability  and  and  control  derivatives.  Based  on 
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Equations  2.3-1  and  2.3-2  the  resulting  dimensional  aerodynamic  force  and 
moment  equations  resolved  into  X,Y  and  Z  axes  components  are: 


Aerodynamic  Force  Components: 


FXa.ro  •  <Xo  ♦  V«*  V  V  V  XTR 

*  XSP'SSP*  "ge^ge1'**. 


Fr«ro*{<VV  V-V*  Vd 

*vr*  v  +  v«s  p)vV4r'V 

FZa.ro  ’  Zo  *  Z«a*  Zi-i*  V  0  *  V«u  * 

Z«e'4e  *  ZTR  *  ZSP4SP*  ZGE'XGE  *  Z«s'4S)'KZg 
Aerodynamic  Moment  Components: 

TXa,n,'(V:CL*LTR>-LGE-’'*LB-»*VP 

+  *T‘r  +  L6#  *  4a+  Lsp'4sp}'xig*  L«r'4r'CL8 


TY  *  (M  +  M  -o+  M. -a+  M  -q  +  M  -6  +  M,  •  5 
aero  '  o  o  o  q  u  u  6  e 

C 

*  mtr  +  MSP'fi$P  +  mge'kge  +  M6s'5S  + 

^aaro  •  {<Vc„  *  ntr>-nge"'  *  Np-p  ♦  Np.r 
*  \  -4a*  »W4f*)-S *  \Sr-\r 


2.3-4 


These  equations  are  conrnonly  referenced  to  either  the  body  axes  or  the 
stability  axes.  If  the  equations  are  referenced  to  the  body  axes,  the 
sln(o)  and  cos(a)  terms  in  the  equations  of  motion  (i.e.  Equations  2.2-1  and 
2.2-2)  are  automatically  set  to  0.0  and  1.0  respectively. 


2.3.2  Nondimensional  Force  and  Moment  Equations 


The  nondimensional  force  and  moment  equations  are  based  on  the  assumption  that 
a  force  (F)  or  moment  (M)  acting  on  any  solid  moving  through  a  fluid  can  be 
expressed  as: 

F  *  CFpl/2  V2A 

M  *  CMpl/2  V2AL 

where:  Cp  *  dimensionless  force  coefficient 

CM  *  dimensionless  moment  coefficient 
p  *  density  of  air 
V  3  airspeed 
A  3  characteristic  area 
L  3  characteristic  length 

The  coefficients  C_  and  Cu  are  assumed  to  be  continuous  functions  of  the  same 

r  n 

variables  (i.e.,  angle  of  attack,  sideslip,  etc.)  used  in  the  dimensional 
case.  For  small  perturbations  about  the  trim  configuration  these  coefficients 
can  also  be  expanded  in  a  Taylor  series  expansion  given  below: 


C 


F 


(4lX;! 


♦  t 


CM 


/Vl 

\wzl'Xt  *  — 


2.3-5 


2.3-6 


The  terms  Cp  and  CM  are  the  aerodynamic  bias  coefficients  for  trim  flight, 
o  o 

The  first  order  derivative  terms  are  nondimensional  stability  and  control 
derivatives.  The  resulting  nondimensional  aerodynamic  force  and  moment 
equations  based  on  Equations  2.3-5  and  2.3-6  resolved  into  components  are: 
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Aerodynamic  Force  Components: 


°e 

♦  Cy  4  Cy  •  icp  *  Cy  •  K»  )*KCy 

*TR  *S P  *  *GE  CSE  8 


"..r.  *  l/2»v2$-  [{(CYs”CCy 


♦  cv  )-Kr  -0 

ytr  cge 


♦  (Cy^-B+  CY  *P  4  C y  - r) •  (b/2V) 

*%'*•*  cvi”)'"VcV1'  c,8r] 


FZmw|  -  l/2pV2S.(C_  4  c7  C-  .-a  ♦  C  .q) 
••ro  z0  za  a  za  Zq 


•(c/2V)  4  C7  '6  +  C7  •  6  4  C7 

4u  u  *  4TR 

e 

4  C7  *  5cp4  ^7  -K,  +  C7  *6^ )*KC7 

ZSP  5P  ZGE  CGE  z5  s  5  ZB 


Aerodynamic  Moment  Components: 


TX^  *  1/2oV2SI>  [((c,  *  4  C,  )  Kr  -6 

aero  U  Lfl  CL  Lyp  CQE 

♦  (Cu.-B+  CL  -p  4  CL  -r) * ( b/2V) 

8  p  r 

♦  C.  •  6  4  C.  •6P.)-KC.  4  C.  '6_’C, 

L4j  a  LfsFSJ  lg  lir  r 


"oro  ’  1««»  SC-(CM#  *  (‘nj**  C^.qJ.tc/ZV) 

♦  Cy  *  U  4  Cy  *  5_4  Gy  4  Cy  *  5c p4  Cy  * 

"u  •  TR  MSP  5P  MGE  CGE 


V  •  l/2ov2sb-  [i%'VV'VB* 

4  (CN..i4  CNp.p  4  CN^.r) i(b/2V)  4  CN  .< 


2.3-7 


2.3-8 
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These  equations  are  commonly  referenced  to  either  the  body  axes  or  the 
stability  axes  reference  systems.  If  the  equations  are  referenced  to  the 
body  axes,  the  sin(a)  and  cos(a)  terms  in  the  equations  of  motion  (i.e. 
Equations  2.2-1  and  2.2-2)  are  automatically  set  to  0.0  and  1.0  respectively. 
The  wing  area  (S)  the  wing  span  (b)  and  the  mean  chord  length  (c)  are  shown  in 
Figure  8.  Figure  9  shows  the  stability  axis  lift  and  drag  forces. 

2.4  Component  Formulation 

The  mathematical  model  for  a  6-D0F  rigid  body  aircraft  are  given  in  Equations 
2.2-1,  2.2-2,  2.2-3,  2.2-4,  2.3-3,  2.3-4,  2.3-7  and  2.3-8.  These  equations 
have  been  implemented  into  the  various  EASY  components  which  are  discussed  in 
the  subsequent  paragraphs. 

2.4.1  Generalized  6D0F  Component  SG 

A  flowchart  for  this  component  is  shown  in  Figure  10.  The  six  degree  of 
freedom  component  SG  determines  the  rates  of  the  aircraft  states  from  the 
rigid  body  equations  of  motion  and  the  kinematic  equations  (Equations  2.2-1 
through  2.2-4).  It  should  be  noted  that  the  equations  2.2-2  were  derived 
assuming  an  x-z  plane  of  symmetry,  whereas  the  equations  contained  in  SG  do 
not  make  this  assumption  and  thus  include  the  contribution  of  the  terms  I 
and  Iyz*  The  parametric  moment  equations  for  component  SG  are: 

TX1  *  TXaero  c°sa-  naero  s1no-  r«P  -  *  Uy-Ij)  qr  *  (pr'P)  !*y 

-  (o2-r2)iy2 

™  *  "«ro  •  lyP  -  Wr2-P2>  -  <Vlx>n>  -  (Pr*P)!xy 

+  (pq-r)Iyz  2.4-1 

TZ1  *  TZaero  cosa  *  TXaerQ  sino*  1 2r  -  I2X(p-qr)  -  (Ix-Iy)pq  -  (P2-d2)Ixy 

-  (pr*-q)Iyz 


17 


Figure  9  Lift  and  Drag  Acting  on  an  Airplane 
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OUTPUTS 


OUTPUTS 


Figure  10  Flowchart  for  Subroutine  SG 
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The  component  receives  as  input  all  body  axis  forces  and  moments  including 
gravity.  Parametric  inputs  include  mass,  moments  of  inertia  and  cross 
products  of  inertia  of  the  aircraft.  Outputs  from  the  component  include  the 
accelerations  and  velocities  (linear  and  angular)  in  the  body  axis  coordinate 
system,  Euler  angles  and  earth  axis  positions.  The  differential  equations  for 
the  twelve  states  representing  the  linear  and  angular  velocities  and  positions 
are  contained  in  this  component  (see  Figure  7). 

2.4.2  Six  DOF  Component  DS 

This  component  is  same  as  the  component  SG  except  that  the  X-Z  plane  is 

assumed  to  be  a  plane  of  symmetry.  Therefore  the  cross  product  of  inertia 

terms  I  and  I  are  zero  and  Equations  2.2-2  apply.  In  addition  the 
xy  y  z 

airplane  positions  in  the  X  and  Y  earth  axis  are  not  evaluated. 

A  flowchart  for  the  component  is  shown  in  Figure  11. 

Lower  DOF  Models: 


The  EASY  program  has  the  capability  of  reducing  the  order  (degrees  of  freedom) 
of  a  system  by  "freezing"  state  variables.  This  is  a  very  useful  technique, 
especially  for  troubleshooting  an  unstable  system.  However,  if  repeated 
analyses  with  reduced  degrees  of  freedom  are  anticipated,  components 
specifically  tailored  to  the  required  degrees  of  freedom  should  be  used.  For 
some  types  of  ACLS  analysis,  for  example  ground  directional  control  studies, 
lateral  stability  analysis,  longitudinal  stability  analysis,  or  drop  test 
correlation,  a  reduced  order  model  can  be  used  with  no  sacrifice  in  accuracy. 
In  such  cases,  there  are  two  significant  advantages  to  using  a  reduced  order 
model,  namely. 

1.  Computational  efficiency  -  a  typical  6  DOF  ACLS  model  can  be  expected 
to  be  about  24th  order.  Using  a  4  DOF  model  where  it  is  appropriate 
reduces  the  order  to  20.  The  saving  in  execution  time  is  typically 
proportional  to  the  square  of  the  model  order,  see  Volume  III, 

Section  2.4. 
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INTEGRATE  LINEAR 
ACCELERATIONS 


LINEAR  ACCELERATIONS 
MOMENTS  (TORQUES) 
MOMENTS  OF  INERTIA 


OUTPUTS 


Figure  II  Flowchart  for  Subroutine  OS 


2.  Clarity  of  output  -  Ease  of  interpretation  of  the  outputs  of  several 
of  the  EASY  linear  analyses  is  considerably  enhanced  if  unnecessary 
states  are  removed  from  the  model . 

2. A. 3  4DQF  Rigid  Body  Dynamics  FD 

The  four  degree  of  freedom  component  consists  of  the  same  types  of  inputs  and 
outputs  as  the  6D0F  component  DS.  Its  degrees  of  freedom  are  the  u  and  v 
components  of  the  linear  velocity  and  the  p  ana  r  components  of  the  angular 
velocity.  This  model  is  applicable  to  ground  directional  control  studies. 

2.4.4  3D0F  Longitudinal  Rigid  Body  Dynamics  TL 

This  is  a  simplified  version  of  the  six  DOF  component  DS  where  only  the 
longitudinal  equations  of  motion  are  retained.  The  degrees  of  freedom  are  the 
u  and  w  components  of  the  linear  velocity  and  the  q  component  of  the  angular 
velocity.  The  model  is  applicable  to  long  period  longitudinal  stability 
analysis  (Phugoid). 

2.4.5  3D0F  Lateral  Rigid  Body  Dyanmics  TD 

This  is  also  a  simplified  version  of  the  six  DOF  component  DS  where  only 
lateral  equations  of  motion  are  retained.  The  degrees  of  freedom  are  the  v 
component  of  the  linear  velocity  and  the  p  and  r  components  of  the  angular 
velocity.  The  model  is  applicable  to  lateral  stability  analysis. 

2.4.6  2D0F  Longitudinal  Rigid  Body  Dynamics  TT 

This  component  is  similar  to  the  component  TL.  The  degrees  of  freedom  are  the 
w  component  of  the  linear  velocity  and  the  q  component  of  the  angular 
velocity.  The  model  is  applicable  to  drop  test  simulations  and  for  analyzing 
inflight  short  period  longitudinal  oscillations. 

2.4.7  Aerodynamic  Variables  Component  VA 

The  model  adds  total  wind  perturbation  onto  the  aircraft  velocity  states, 
and  computes  aero  variables  such  as  angle-of-attack,  sideslip,  dynamic 
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pressure,  and  mach  number.  The  cos(o)  and  sin(a)  for  the  body  to  stability 
axis  transformation  is  calculated  in  VA.  A  switch  is  used  to  assure 
compatabi 1 ity  with  various  DOF  components.  This  concept  is  illustrated  in 
Figure  12. 

2.4.8  Longitudinal  Forces  and  Moments  Component  OL 

Inputs  to  this  component  include  external  longitudinal  forces  and  torques 
(e.c.  due  to  trunk,  engine)  and  the  longitudinal  aero  stability  derivatives 
used  to  calculate  the  longitudinal  aerodynamic  forces. and  moments.  The 
external  forces  and  moments  are  combined  with  the  calculated  aerodynamic 
forces  and  moments  in  body  axes  and  output  as  the  longitudinal  force  and 
moment  vectors  about  the  aircraft  c.g.  The  body  axis  linear  accelerations  are 
determined  along  with  the  moments  which  are  passed  to  the  equations  of  motion 
module  DS  or  SG. 

The  aero  derivatives  can  be  specified  in  either  stability  axes  or  in  body 
axes.  For  more  details  see  Section  5.1  of  the  Users  Manual.  A  flowchart  for 
the  component  OL  is  shown  in  Figure  13. 

2.4-9  Lateral  Forces  and  Moments  Component  DL 

This  component  is  identical  in  design  to  component  OL  and  computes  the  lateral 
set  of  linear  accelerations  and  moments.  A  flowchart  for  the  component  is 
shown  in  Figure  14. 

2.5  Wind  Components 

2.5.1  Gust  Wind  Model  (GW) 

The  gust  wind  model  (GW)  is  a  continuous  random  turbulence  model  of  the  Dryden 
form.  Background  information  and  a  description  of  the  model  are  given  in  the 
next  two  subsections. 
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Figure  12  Use  of  Switches  for  Various  DOF  Models 
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RMCIS  AND  MHCNTS.  ACM  VA  A! ASICS  ANO  CONTAOL  SUAFACES 
KFLCCTIONS 


INPUTS 


2.S.1.1  Background 


The  gust  wind  model  was  taken  from  MIL-F-8785B,  Section  3.7.  In  general,  the 
gust  wind  models  are  classified  into  continuous  and  discrete  types.  This 
military  specification  presents  two  continuous  gust  wind  mathematical  models, 
the  Von  Karman  form  and  the  Dryden  form.  The  Dryden  form  was  implemented  in 
GW  because  the  Von  Karman  form  was  designed  to  be  compatible  with  aircraft 
Structural  analysis  which  is  of  no  interest  in  a  rigid  body  model. 

The  Dryden  form  of  the  gust  wind  spectra  for  linear  and  angular  turbulent  body 
axes  velocity  components  expressed  in  transfer  function  form  are: 

Linear  Velocity  Components: 
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Angular  Velocity  Components: 
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The  scales  and  intensities  for  clear  air  turbulence  and  thunderstorm 
turbulence  given  in  MIL-F-8785B  are  reproduced  here  to  familiarize  the  reader 
with  typical  values  for  these  variables.  The  clear  air  turbulence  given  in 
this  MIL  Standard  defines  a  root-mean-square  intensity  for  ow  as  a  function  of 
altitude,  see  Figure  15.  The  clear  air  intensities  and  cy  may  be  obtained 
from  the  following  relationship: 


2.5-3 


The  clear  air  turbulence  scales  needed  to  evaluate  the  above  equations  are 
divided  into  the  two  altitude  (h)  ranges  given  below: 
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For  h  above  1750  feet: 

L  ■  L  ■  L  3  1750  feet 
u  v  w 

For  h  below  1750  feet: 

L  ■  altitude 

W 

L  •  L  «  145  h1/3  feet 
u  v 

The  scales  and  intensities  for  a  typical  thunderstorm  given  in  this  MTL 
Standard  are  listed  below: 

a  *  o  3  o  3  21  ft/sec 
u  v  w 

L  ■  L  3  L  *  1750  feet 
u  v  w 

These  are  two  examples  of  typical  intensities  and  scales  designed  for  the 
Dryden  form  of  continuous  turbulence  model.  The  user  should  be  aware  that 
there  are  other  statistical  curves  which  have  also  been  developed  for  use  with 
this  model. 

2.5. 1.2  Model  Description 

The  gust  wind  model  implements  the  Dryden  form  for  the  linear  and  angular 
transfer  function  velocity  components  give  in  Equations  2.5-1  and  2.5-2.  Also 
the  components  allows  Ug,  Vg,  wg  and  pg  to  be  modified  by  random  noise  samples 
generated  by  either  subroutines  RA  and  RN.  To  implement  GW  the  user  must 
specify  the  characteristic  intensities  o^,  oy  and  ow  and  the  scales  Lu,  Ly  and 

Lw* 

A  block  diagram  of  a  complete  wind  turbulence  model  incorporating  subroutines 
RN  and  RA  along  with  the  gust  wind  (GW)  and  shear  wind  (WS)  components  are 
shown  in  Figure  16. 

2.5.2  Shear  Wind  Model 

The  model  simulates  the  boundary  layer  reduction  of  wind  velocity  as  the 
airplane  approaches  the  ground.  The  velocities  output  from  the  shear  wind 
model  are  combined  with  any  gust  or  turbulence  components  to  determine  net 
air-stream  velocity  components,  see  Figure  16. 
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Figure  16  Wind  Turbulence  Model  Schematic 


The  wlndshear  velocity  is  calculated  by  multiplying  the  wind  magnitude  quoted 
at  control  tower  altitude  of  50  feet  by  a  non-linear,  dimensionless  shear 
factor,  MF,  which  reflects  the  change  in  wind  with  altitude.  Thus 

MF  •  f  (Altitude)  2.5-5 

MF  is  a  user  defined  tabular  input.  Two  sample  profiles  are  shown  In  Figure 
17. 

The  wind  magnitude  vector,  WK,  (ft/sec)  is  assumed  parallel  to  the  ground 
plane  at  an  angle  WAN  from  the  centerline  of  the  runway.  The  two  windshear 
components  along  and  perpendicular  to  the  runway  are: 

WKN  -  -WK*MF*CQS  (WAN) 

WKE  *  WK*MF*$IN  (WAN)  2.5-6 


NORTH 
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MEAN  WIND  VELOCITY 


Figure  17  Wind  Shear  Prof 11 


These  two  runway  components  are  then  transformed  into  X,Y,Z  body  axis 
components  UWS,  VWS  and  WWS  respectively. 

UWS  =■  WKN*C0S  (Pitch) 

VWS  =■  WKE  2.5-7 

WWS  =  -WKN*SIN  (Pitch) 

2.6  Engine  Model s 

Two  engine  components  were  developed.  Both  calculate  engine  thrust  and  engine 
bleed  and  fan  air  pressures  and  temperatures.  The  two  modules  differ  only  in 
the  way  the  engine  dynamic  effects  are  simulated  as  explained  below. 

2.6.1  First  Order  Lag  Model,  ES 

The  model  ES  is  dynamically  simple,  with  all  dynamic  effects  (e.g.  engine  spin 
up  and  spin  down)  accounted  for  by  a  first  order  lag  for  the  thrust 
calculation. 


THRUST 


A  flowchart  for  the  component  is  shown  in  Figure  18.  The  thrust  output  TH  is 
a  state  variable  an  is  computed  by  the  following  equation: 

TH  =  (THR-TH)/TC0 

where: 

THR  ■*  Required  Thrust  Level,  lbs. 

TCO  ■  Engine  Time  Constant,  sec 

TH  =  Time  rate  of  change  of  TH,  Ib/sec 

TH  =  Thrust  at  current  time,  lbs 

Any  external  force  affecting  the  thrust  (e.g.  engine  powered  thrusters)  is 
then  added  to  TH  and  the  net  thrust  is  transformed  into  vehicle  body  axes 
forces  and  moments. 

FX  *  TH*GAX 
FI  =  TH*GAZ 
TY  *  Z0*FX-X0*FZ 

where: 

FX,  FZ  =  X,  Z  body  axis  forces,  lbs 

TY  =  Y  body  axis  torque  (pitching  moment),  ft-lb 

GAX,GAZ  =  X,  Z  body  axis  direction  cosines 

XO,  ZO  =  X,  Z  engine  moment  arms,  ft. 

These  forces  and  moments  are  then  supplied  to  the  longitudinal  aerodynamic 
force  and  moment  component  OL. 

The  remaining  portions  of  the  model,  below  the  dotted  line  in  Figure  18, 
calculate  the  bleed  and/or  fan  air  properties  if  requested.  The  bleed  air 
properties  (P,T,W)  would  be  needed  for  example  tf  bleed  air  powered  the  an  air 
cushion  system  (trunk).  Similarly  if  the  engine  was  a  turbofan  type,  the  fan 
stage  properties  may  be  of  interest.  The  user  option  is  indicated  by  flag 
indicators  (user  inputs).  The  following  equations  are  used: 


Figure  18  Flowchart  for  Subroutine  ES 


Engine  Inlet  Ram  Rise: 


Total  pressure  (PT)  and  temperature  (TT)  are  calculated  using  standard 
compressible  flow  equations: 


where : 


pT  *  Wd*-^)3'5 


M 

p 

a  mb 
Tamb 


*  Mach  No. 

3  Ambient  Pressure,  psia 
3  Ambient  temp.,  °R 


Engine  Speed: 


SPO  =  f  (M,  TH)  See  Figure  19  for  typical  user  input. 

ENC  3  SPD  ^519/TT 

where: 

ENC  *  corrected  engine  speed 


Bleed  Air  Properties: 

Compression  Temperature  rise  3  DT 

DT  *  f  (ENC)  See  Figure  20  for  typical  user  input. 

Tport  ’  TT*(U0T>  ’  TPU 
Compression  Pressure  Ratio  3  CPR 

CPR  3  f  (ENC)  See  Figure  21  for  typical  user  input. 

Bleed  Port  Upstream  Pressure  3  PPU  3  PT*CPR 
Port  Pressure  Ratio  3  PRAT 
PRAT  3  PPU/P2 

where  P2  3  port  downstream  pressure,  psia 
Bleed  Air  corrected  flow  rate  3  WCR 

WCR  3  f(PRAT)  See  Figure  22  for  typical  user  input. 

“Bleed  "  W2‘  lb/min 
W2  -  WCR*PPU/  TPU 

T2  3  TPU 
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RELATIVE  SPEED 


Figure  21  Engine  Compressor  Pressure  Ratio 
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Figure  22  Engine  Port  Pressure  Drop  Characteristics 
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The  fan  air  properties  are  calculated  using  similar  equations  for  pressure  and 
temperature  rise.  The  bleed  and  fan  air  property  calculations  are  identical 
in  both  engine  models. 

2.6.2  Second  Order  Transfer  Function  Engine  Model,  EC 

This  more  complex  model  (EC)  is  shown  in  Figure  23  in  block  diagram  form.  The 
bleed  and  fan  air  properties  calculations  are  identical  to  those  of  component 
ES  and  therefore  are  not  discussed  here.  The  engine  pressure  ratio  (EPRI)  is 
computed  as  a  function  of  the  power  iever  angle.  Engine  transient  conditions 
are  accounted  for  by  introducing  a  first  order  lag  for  engine  deceleration  and 
a  second  order  lag  for  engine  acceleration.  Net  thrust  (forward  or  reverse) 
is  computed  as  a  function  of  mach  number  and  the  actual  engine  pressure  ratio, 
EPRO.  The  second  order  function  requires  user  inputs  for  engine  spinup 
natural  frequency  versus  f(EPRO)  and  damping  ratio  constants.  The  following 
equations  represent  the  engine  model  described  in  block  diagram  form  in 
Figures  23  and  24. 


JUATIONS: 


EPRI  *  f(dT)  Figure  25a 

where  0-j.  =  power  level  angle 
AE  *  EPRI -EPRO 
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Figure  23  Engine  Model  EC 
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Figure  25  Typlcel  Jet  Engine  Char«cttr1«t1ei 


Solving  for  EPRO  and  EPRI  in  terms  of  w^  gives: 

EPRO  a  _ 

EPRI  s2  +  2^wns  +  w2 

Engine  Forward  Thrust: 

Tp  -  f  (EPRF)  since  EPRF=EPR0 
Engine  Reverse  Thrust: 


Tr  «  -f  (EP RR) 


where: 


EPRR  = 


EPRO 

T+Ts 


Figure  25c 


Figure  25d 


2.7  Thruster  Models 

Due  to  dynamic  instability  experienced  during  ground  handling  on  the  ACIS 
equipped  aircraft,  auxiliary  control  thrusters  were  proposed  to  alleviate 
this  problem  (see  AFFDL -TR -77-21 ,  Final  Report  on  Jindivik  ACLS  Integration 
Contract  F33615-75-C-3088) .  Individual  thrusters  were  recomnended  for  each  of 
the  yaw,  pitch  and  roll  modes.  The  present  contract  required  that  specific 
provision  be  made  to  simulate  the  effects  of  auxilliary  yaw,  pitch  and  roll 
thrusters. 

2.7.1  Yaw  Control  Thruster,  YC 

The  purpose  of  component  YC  is  to  simulate  the  effects  of  the  Y  body  axis 
thrust  produced  by  an  auxilliary  yaw  thruster.  The  yaw  control  thruster 
component  will  model  either  an  aircraft  engine  powered  thruster  or  a 
self-powered  thruster. 


i 
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A  block  diagram  flow  chart  of  the  model  is  shown  in  Figure  26.  Various 
input/output  variables  are  listed  below: 


*****  Outputs  ***** 

FX  Engine  Thrust  Reduction,  lbs  (0  =  if  self  powered  thruster  is  modeled) 

FY  Vectored  Thrust-Side  Force,  lbs 

TX  Roll  Moment  Due  to  Thruster,  ft- lb 

TZ  Yaw  Moment  Due  to  Thruster,  ft-lb 

*****  inputs  ***** 

ED  Engine  Dependence  Indicator(ED=1.0, YES  ;ED*0.0,N0) 

TM  Thruster  Maximum  Force  for  Engine  Independent  System,  -lbs 

ST  Slope  for  Maximum  Available  Side  Thrust  as  Function  of 

Engine  Thrust 

SR  Slope  of  Engine  Thrust  Reduction  as  Function  of  Vectored 

Thrust 

Cl  Saturation  Function  Slope 

C2  Saturation  Slope 

SIG  Aircraft  Control  System  Signal  to  Thruster  (e.g.  Yaw 

Angle) 

GA  First  Order  Lag  Gain 

TC  First  Order  Lag  Time  Constant,  seconds 

TH  Engine  Thrust,  lbs 

XA  Thruster  Yaw  Moment  Arm,  ft 

ZA  Thruster  Roll  Moment  Arm,  ft 


The  equations  used  to  model  the  yaw  control  thruster  are  shown  in  the  flow 
chart. 


2.7.2  Pitch  Control  Thruster  PT 


The  pitch  control  thruster  produces  a  body  axis  thrust  which  augments  pitch 
stability  of  an  ACLS  equipped  vehicle.  The  component  is  identical  in  design  to 
YC  and  a  separate  flowchart  for  the  component  PT  is  unneccessary.  The 
input-output  variable  list  is  given  below: 


*****  Outputs  ***** 

FX  Engine  Thrust  Reduction 

FZ  Vectored  Thrust-Vertical  Force 

TY  Pitch  Moment  Due  to  Thruster 

*****  Inputs  ***** 

ED  Engine  Dependence  Indicator  (ED*1.0,YES;ED*0.0,N0) 

TM  Thruster  Maximum  Force  for  Engine  Independent  System 

ST  Slope  for  Vectored  Thrust  as  Function  of  Engine  Thrust 

SR  Slope  of  Engine  Thrust  Reduction  as  Function  of 

Vectored  Thrust 

Cl  Saturation  Function  Slope 
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YAW  CONTROL  THRUSTER  (YC) 


FI gurf  26  Flowchart  for  Sub rout 1 no  YC 


C2  Saturation  Slope 

SIG  Aircraft  Control  System  Signal  to  Thruster 

GA  First  Order  Lag  Gain 

TC  First  Order  Lag  Time  Constant 

TH  Engine  Thrust 

XA  Thruster  Pitch  Moment  Arm 


2.7.3  Roll  Control  Thruster  RT 


The  roll  control  thruster  produces  a  Z  body  axis  thrust  which  augments  roll 
stability  of  an  air  cushioned  vehicle.  The  roll  control  thruster  is  also 
identical  in  design  to  the  other  two  thrusters,  YC  and  PT.  The  input-output 
variable  list  is  given  below: 


*****  Outputs  ***** 

FX  Engine  Thrust  Reduction 

FZ  Vectored  Thrust-Vertical  Force 

TX  Roll  Moment  Due  to  Thruster  Not  on  X-axis 

*****  inputs  ***** 

ED  Engine  Dependence  Indicator(ED*1.0,YES=0.0,N0) 

TM  Thruster  Maximum  Force  for  Engine  Independent  System 

ST  Slope  for  Vectored  Thrust  as  Function  of  Engine  Thrust 

SR  Slope  of  Engine  Thrust  Reduction  as  Function  of  Vectored 

Thrust 

Cl  Saturation  Function  Slope 

C2  Saturation  Slope 

SIG  Aircraft  Control  System  Signal  to  Thruster 

GA  First  Order  Lag  Gain 

TC  First  Order  Lag  Time  Constant 

TH  Engine  Thrust 

YA  Thruster  Roll  Moment  Arm 


SECTION  III 


AIRFLOW  COMPONENTS 


3.1  Ducting  Components 

Components  included  under  the  general  heading  of  ducting  components  include 
ducts,  splits,  merges,  and  valves.  Equations  used  are  for  the  simulation  of 
the  dynamic  flow  of  gases.  Bends,  expansions  and  contractions  are  not 
considered  as  separate  components,  since  subdividing  a  length  of  ducting  into 
many  small  lengths  gives  rise  to  very  high  frequency  terms  in  the  dynamic 
simulation  which  have  no  practical  significance.  The  pressure  drop  effects  of 
bends  etc.  can  however  be  included  in  the  so  called  k  factor  input. 


3.1.1  to  3.1.4,  DU,  DV,  FS,  MG  &  FR  Components 


Under  a  previous  Air  Force  contract  for  Environmental  Control  System  (ECS) 
several  airflow  components  were  developed  and  fully  documented  (see  Reference 
4).  Five  of  these  components,  namely,  duct  ( DE ) ,  split  (SP),  merge  (ME), 
valve  in  a  duct  (VD),  and  compressor  with  surge  analysis  (CR)  were  transferred 
from  the  existing  ECS  1 ibrary  to  the  EASY-ACLS  library.  These  components 
contained  specific  heat  and  liquid  water  contents  as  parameters  which  were 
removed  since  they  are  of  no  particular  interest  in  an  ACLS  airflow  system. 

The  modified  components  were  renamed  as  follows: 


COMPONENT 
DUCT  (DE) 

FLOW  SPLIT  (SP) 

FLOW  MERGE  (ME) 
COMPRESSOR /FAN  WITH 
SURGE  ANALYSIS  (CR) 
VALVE  IN  A  DUCT  (VD) 


ACLS  LIBRARY  NAME 
DU 
FS 
MG 


FR 

DV 


SO 


Since  all  necessary  documentation  of  these  components  are  available  in  tne 
Reference  4,  Volume  I  no  additional  documentation  is  provided  here. 


3.2  Fan  and  Ejector  Components 
3.2.1  Ejector,  EJ 

The  purpose  of  the  component  EJ  is  to  model  an  air  ejector  with 
converging-diverging  nozzle  and  subsonic  or  choked  flow  conditions.  The 
subroutine  uses  a  two  dimensional  input  table  of  flow  ratio  (total /primary)  as 
a  function  of  the  two  pressure  ratios  (total /secondary)  and  (primary/ 
secondary).  For  choked  throat  flow  upstream  pressure  is  computed  to  match 
flow  and  for  subsonic  flow  the  exit  pressure  equals  secondary  supply  pressure. 


P2-W2 


\ 

P  W 
3  3 
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The  input-output  and  state  variable  names  used  by  component  EJ  are  defined 
below: 


*****  Inputs  ***** 

T1  -Inlet  Temperature  Primary  Source,  Oeg  Rankine 

W1  -Inlet  Flow  Rate  Primary  Source,  lbs/min  (Port  No.  1) 

T2  -Inlet  Temperature  Secondary  Source,  Deg  Rankine 

P2  -Inlet  Pressure  Secondary  Air  Source,  psia  (Port  No.  2) 

P3  -Outlet  Pressure,  psia  (Port  No.  3) 

ANT  -Nozzle  Throat  Area,  sq.  ft. 

ANE  -Nozzle  Exit  Area,  sq.  ft. 

AK  -Convergent-Divergent  Nozzle  Diffuser  Loss  Factor 

-(For  Convergent  Nozzle,  Input  AK=0.,  ANE=ANT) 

*****  Table  ***** 

TAB  -Flow  Ratio(Total /Primary)  as  a  Function  of  the  TV/C 

Pressure  Ratios(Total /Secondary  and  Primary/Secondary) 
Two  Dimensional  Table 
*****  Outputs  ***** 

T3  -Outlet  Temperature,  Deg  Rankine  (Port  No.  3) 

W3  -Total  Outlet  Flow,  lb/min  (Port  No.  3) 

PI  -Inlet  Pressure  Primary  Air  Source,  psia  Port 

P1D0T  -Inlet  Pressure  Derivative,  psia/sec  No.  1 

We  can  now  write  the  following  equations  defining  the  ejector  model. 


If  W1  is  greater  than  WCHO,  choked  throat  flow  is  present. 

WCHO  =  31.9*ANT*P1/SQRT(T1) 

Choked  Flow  -  Calculate  Upstream  Pressure  to  Match  Flows 
P1CAL  =  W1*SQRT(T1)/(31.9*ANT) 

P1D0T  =  (P1CAL-P1)/.01 

Flow  Is  Not  Choked  -  Calculate  Dynamic  Head  at  Throat 
AM  =  AMACH (PI ,T1 ,ANT,W1 ,0. ) 

GAMMA=1 .4 

PTS  =  Pl/(1 .+(GAMMA-1 ) *AM* AM/2. **( GAMMA/ (GAMMA-1 . ) ) 

ALOSS  =  (Pl-PTS )*AK 

Calculate  Total  Pressure  at  Nozzle  Exit 
PE  *  PI -ALOSS 

AME  =  AMACH (PE ,T1 ,ANE ,W1 ,0. ) 

PESCAL  =  PE/( 1 .+(GAMMA-1 . )*AME*AME/2. )**( GAMMA /( GAMMA- 1. ) ) 
P1D0T=(P2-PESCAL)/.01 

Knowing  P3/P2  and  P1/P2,  WRAT  =  ( W>/T^ 3/ ( wVt) ^  can  be  determined  from  Table  TAB. 
Now,  T3  may  be  calculated  using  the  energy  balance  equations: 


W3-T3  *  Wl-Tl  +  W2-T2 


and  by  assuming  a  steady  state  exists  which  implies: 


W3  -  W1  +  W2 
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The  theory  for  fans  and  compressor  models  was  developed  in  Reference  4,  Volume 
I,  Section  3.2.2. 

3.2.2. 1  Centrifugal  Fan,  FR 

The  module  FR  is  an  adapted  version  of  the  ECS  module  CR  (See  reference  4, 
Section  3.2.2).  The  model  calculates  the  component  performance  using  steady 
state  maps.  Inlet  pressure  is  the  one  state  variable  for  the  model.  The 
model  is  capable  of  reverse  or  surge  flow  analysis  and  can  accomodate  pressure 
ratios  of  less  that  1. 

3. 2.2.2  Axial  Fan  (with  Hysteresis),  FH 

The  purpose  of  the  component  FH  is  to  model  an  externally  driven  axial  fan 
with  transition  time  constant  between  stall  and  recovery.  It  is  assumed  that 
the  transition  from  normal  to  stalled  operation  occurs  when  the  pressure  ratio 
exceeds  the  stall  pressure  ratio  PRS  (user  supplied)  and  the  transition  from 
stalled  to  normal  operation  occurs  when  pressure  ratio  falls  below  a  reverse 
pressure  ratio  PRR  (user  supplied). 
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The  inputs  (including  tables),  outputs  and  state  variables  are  listed  below. 
The  fan  flow  rate  and  the  inlet  pressure  are  two  states. 

*****  Tables  *  *  *  *  * 

CF  -  Fan  Flow  Rate(Forward)  as  a  Function  of  Fan 
Pressure  Ratio,  One  Dimensional  Table 
CR  -  Fan  Flow  Rate(Reverse)  as  a  Function  of  Fan 
Pressure  Ratio,  One  dimensional  Table 
*****  Outputs  ***** 

WC  -  Fan  Flow  Rate,  Ib/min 
WCDOT  -  Fan  Flow  Rate  Derivative,  Ib/min/sec 
PI  -  Inlet  Pressure  (Port  No.  1),  psia 
P1D0T  -  Inlet  Pressure  Dertivative,  psia/sec 
T2  -  Outlet  Temperature  (Port  No.  2),  Deg  R 
*****  Inputs  ***** 

P2  -  Outlet  Pressure  (Port  No.  2),  psia 
W1  -  Inlet  Flowrate  (Port  No.  1),  lb/mi n 

T1  -  Inlet  Temperature  (Port  No.  1),  Deg  R 

PRR  -  Pressure  Ratio  Below  Which  Transition  From 
Stalled  to  Normal  Operation  Occurs 
PRS  -  Pressure  Ratio  Above  Which  Transition  From 
Normal  to  Stalled  Operation  Occurs 
TC  -  Fan  Time  Constant,  sec 
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The  main  equations  used  are: 


T2  3  T1 

PI CAL  *  W1*SQRT(T1)/WC 
PI DOT  *  (P1CAL-P1 )/.01 
PR  -  P2/P1 
WCF  *  CF(PR) 

WCR  *  CR(PR) 

WCCAL  3  WCF  if  WOWCF  &  PR<PRR 

WCCAL  3  WCR  if  WC<WCR  4  PR>PRS 

WCDOT  3  (WCCAL -WC)/TC 


3. 2. 2. 3  Axial  Turbofan,  FT 


The  purpose  of  fan  component  FT  is  to  simulate  a  hub  or  tip  driven  axial 
turbofan  as  used  on  the  Jindivik  ACLS  vehicle.  The  model  uses  input  tables 
defining  steady  state  characteristics  of  turbine  (Figure  27)  and  fan  (Figure 
28).  Drive  air  turbine  inlet  pressure  is  a  state  variable.  The  input,  output 
list  is  as  follows: 


*****  Table  ***** 

WC  -Table  of  Corrected  Turbine  Flow  as  a  Function  of 
-Orive(Bleed  Air)  to  Cushion/Trunk  Pressure  Ratio, 

-One  Dimensional  Table 

TOT  -Table  of  Total  Flow  From  Turbofan  as  a  Function  of 
-Cushion/Trunk  Pressure  (PSFG)  and  Drive  Pressure  (PSIA) 

-Two  Dimensional  Table 
*****  Outputs  ***** 

T3  -Temperature  of  Fan  Air  Exit,  Deg  Rankine 

W3  -Total  Flow  From  Turbofan  to  Cushion/Trunk,  lb/min  (Port  No.  3) 

PI  -Drive/Bleed  Air  Pressure,  PSIA  (Port  No.  1) 

P1D0T  -Derivative  of  PI,  PSIA/sec 
*****  Inputs  ***** 

T1  -Orive/Bleed  Air  Temperature,  Deg  Rankine 
W1  -Orive/Bleed  Air  Flow  Rate,  lb/min  (Port  No.  1) 

T2  -Ambient  Air  Temperature,  Deg  Rankine 
P2  -Ambient  Air  Pressure,  PSIA  (Port  No.  2) 

P3  -Pressure  of  Fan  Air  Exit,  PSIA  (Port  No.  3) 

VOL  -Internal  Volume,  cu.  ft. 

FC  -Frequency  Control  on  PI  (FC.GE.l.) 

-A  Value  of  FC  Greater  Than  1.  Decreases 
-Frequency  Response  of  PI  Correspondingly 


TURBINE  CORRECTED  MASS  FLOW  <*»V4-  .  LB/SEC 


SOURCE  :  AFFDL-TR-77-21  INTEGRATION  OF  ACLS  TECHNOLOGY 
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Figure  27  ACLS  Turbofan  Turbine  Flow  Characteristics 
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TOTAL  CORRECTED  NASS  FLOM  (LBN/SEC) 

[dependent  variable] 


CORRECTED  STATIC  PRESSURE  (PSFS) 
[PRIMARY  INDEPENDENT  VARIABLE] 


Figure  28  ACLS  Fan  Performance  Map 
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The  following  equations  define  the  fan  model. 


PRAT  «  P1/P3  if  >1. 

WCOR  *  f(PRAT) ;  Table  WC,  Figure  27 

W1CAL  =  50*WC0R*(SQRT(519)/14.7 )*P1/SQRT (Tl) 

PI DOT  =  R*T1*(W1-W1CAL)/(60.*144.*FC*V0L) 

PSF  *  (P3-P2)*144 

W3  *  60*f(PSF ,P1 )  if  >W1CAL  Table  TOT,  Figure  23 

W2  *  W3-W1CAL 

T3  *  (W1*T1+W2*T2)/W3  if  >400. 


The  concept  of  parameter  FC  is  based  on  the  ducting  component  theory,  see 
Reference  4,  Volume  I. 


3.2.2. 4  Inlet  Fan 

This  fan  component  is  designed  to  be  positioned  at  the  inlet  of  an  air  flow 

network.  The  user  must  input  the  ambient  temperature  (Tamb),  ambient  pressure 

(Pamb)»  maximum  ram  temperature  (Tram) ,  maximum  ram  pressure  (pram)  and  ram 

efficiency  (n  )•  The  inlet  ram  effects  determine  the  fans  inlet  temperature 
'  ram 

and  pressure,  T  and  P.  given  below: 
in  in  J 

^"in  3  ^ram  “  ^amb^’  ^  ram  +  ^amb 

P.  s  (P  •  P  \  •  n  4*  P 

in  '  ram  amb'  'ram  amb 


The  fan  flow  rate  (Wv&/<5)is  determined  from  a  fan  performance  map  input  as 
functions  of  the  pressure  ratio  (pquj/Pjn)  and  the  ^an  speed.  The  fan  output 
temperature  (Tgyj)  is  determined  by  the  following  equations: 


T 


out 


»  T.  +  AT 
in 
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where: 


Y-l 


R  *  gas  constant 

Cp  ■  specific  heat  at  constant  pressure 

Hfan  ■  fan  efficiency 
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SECTION  IV 


BOEING  INELASTIC  TRUNK  ANO  CUSHION  COMPONENT 


4.1  General  Model  Description 

4.1.1  Trunk  Geometry 

The  model  represents  the  geometry  and  pneumatic  characteristics  of  an 
inelastic  trunk  enclosing  a  cushion  of  air  (Figure  29). 

The  trunk  is  described  by  a  finite  number  of  elements  which  need  not  be 
uniform  in  cross  section,  width,  or  location  on  the  vehicle  (Figure  30). 

The  trunk  is  composed  of  three  parts:  side  elements,  fore  elements  and  aft 
elements.  The  side  elements  have  unconstrained  outward  movement  subject  to 
ground  contact  or  changes  in  the  cushion  trunk  pressure  ratio  and  are 
represented  with  a  membrane  model  as  shown  in  Figure  31.  The  fore  and  aft 
elements  experience  little  or  no  outward  motion  due  to  constraining  peripheral 
stresses  and  are  represented  by  a  frozen  model  or  a  constrained  membrane  model 
(user  option)  as  shown  in  Figure  32. 

The  element  locations  are  specified  in  the  trunk  coordinate  system  by  the 
position  of  the  trunk  inboard  attach  point  (Figure  33).  The  trunk  axis  is 
referenced  to  the  vehicle  body  axis  as  shown  in  Figure  29.  Thus  the  entire 
trunk  assembly  can  be  repositioned  on  the  vehicle  by  changing  the  value  of  two 
parameters,  BST  and  WLT. 

4.1.2  Pneumatic  System 

A  schematic  of  the  trunk-cushion  pneumatic  system  is  shown  in  Figure  34. 
Compressible  flow  equations  are  used  to  describe  the  state  of  the  fluid.  The 
remaining  air  supply  system  not  shown,  such  as  bleed  air,  ejectors,  fans,  etc. 
are  modeled  with  other  EASY  components.  Air  can  be  supplied  to  the  trunk 
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AIR  SUPPLY  TO  TRUNK  AIR  SUPPLY  TO  CUSHION 


-  TRUNK  TO  CUSHION  FLOW  THRU  TRUNK  PERFORATIONS 

-  TRUNK  TO  ATMOSPHERE  FLOW  THRU  TRUNK  PERFORATIONS 
Wra  -  TRUNK  TO  ATMOSPHERE  FLOW  THRU  RELIEF  VALVE 

-  CUSHION  TO  ATMOSPHERE  FLOW  THRU  TRUNK-GROUND  GAP 

Figure  34  Trunk -Cushion  PneumaMc  System  Schematic 
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and/or  cushion.  Air  outflow  occurs  through  the  trunk  perforations,  the 
trunk-ground  gap  and  the  trunk  relief  valve. 

4.1.3  Component  Structure 

A  simplified  flow  chart  of  the  trunk  component  is  shown  in  Figure  35  which 
describes  the  basic  function  of  the  model.  Input  variables  required  are 
velocities  and  displacements  from  the  vehicle  rigid  body  component  and 
flow-rate  and  temperature  from  the  air  supply  components.  During  the  initial 
pass  through  the  component  (i.e.  when  T=0)  parameter  data  tables  are 
calculated  from  a  set  of  implicit  equations  which  represent  the  trunk  section 
properties.  Subsequent  evaluation  of  the  element  section  properties  can  then 
be  made  economically  by  interpolating  these  tabular  data.  A  "DO"  loop  is 
entered  within  the  code  which  evaluates  the  individual  element  properties  and 
calculates  their  accumulative  effect  as  a  trunk-cushion  system.  Note  that  the 
terrain  model  is  called  from  the  trunk -cushion  model  and  is  not  an  independent 
EASY  component. 

Vehicle  forces  and  moments  are  calculated  as  well  as  pressure  and  volume  rates 
of  change  and  flow  rates  which  are  outputted  to  their  corresponding  components. 

4.2  Discussion  of  Model  Idealization 

A  detailed  flow  chart  of  the  trunk-cushion  component  is  shown  in  Figure  36  and 
corresponds  to  the  simplified  chart  of  the  model  shown  in  Figure  35.  The 
detailed  description  of  the  approach,  assumptions  and  equations  of  the  model 
subcomponents  will  be  presented  in  the  same  order  as  they  are  shown  on  the 
flow  chart. 

4.2.1  Trunk  Section  Properties 

Each  element  of  the  trunk  is  represented  by  a  cylindrical  inelastic  membrane. 
The  set  of  equations  describing  its  shape  are  implicit  and  require  an 
Iterative  solution.  Solving  these  equations  in  this  manner  for  each  call  from 
the  EASY  integrator  would  be  very  time  consuming  and  impractical.  The 
approach  used  for  this  model  is  to  initally  solve  the  set  of  equations  (i.e., 
at  T-0)  for  n  values  of  Pr  and  ZQ  and  store  the  parameter  data  in  arrays  (see 
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Figure  35  Trunk  and  Cushion  Component  (Concluded) 


} 


Figure  36  Trunk  Cushion  Component  "IK"  Flowchart  (Concluded) 
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Figure  37).  The  section  properties  are  then  evaluated  from  the  arrays  using 
table  look-up  routines-  during  the  model  execution  (see  Figures  38  and  39). 
The  section  parameter  data  is  generated  by  subroutine  IC  as  shown  in  Figure 
37.  The  model  will  generate  and  store  up  to  six  sets  of  parameter  data  each 
representing  an  element  of  size  and  shape  specified  by  the  user.  Each 
parameter  set  can  be  generated  as  a  membrane  model  or  a  frozen  model . 

4.2. 1.1  Side  Element  (Membrane  Model) 

The  free-body  diagrams  of  the  trunk  element  membrane  and  the  trunk  geometry 
for  loaded  and  free  shape  are  shown  in  Figure  40.  The  following  assumptions 
were  observed  in  the  derivation  of  the  governing  equations. 

1.  The  trunk  element  behaves  as  an  inelastic  membrane. 

2.  Trunk  membrane  tension  is  constant  along  its  minor 
periphery  (i.e.,  no  contribution  from  ground  tractions, 
etc.) 

3.  Pressure  change  from  Pc  to  Pfl  occurs  at  the  venturi 
throat . 

4.  The  trunk  membrane  is  tangent  at  venturi  throat. 

5.  Perforation  nozzle  reactions  are  negligible. 

6.  The  pressures  on  both  sides  of  membrane  element  in 
contact  with  the  ground  are  equal  to  P^. 

7.  The  portion  of  membrane  in  contact  with  the  ground  is  a 
straight  line. 

8.  The  trunk  is  tangent  at  the  edge  of  ground  contact. 

9.  The  pressure  change  from  Pc  to  Pt  and  from  Pt  to  Pa 
occurs  instantaneously  at  the  edge  of  the  flattened 
area. 

The  following  derivation  is  for  the  free  trunk  shape  shown  in  Figure  40a. 
From  the  membrane  free  body  diagram 

*V<4V  2  <Pt-Pc)  r2  sl4  *2 
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Figure  37  Subroutine  "IC"  Flowchart  (Continued) 


Flour*  37  Subroutine  "IC"  Flcwcturt  (Continued) 


C  PURPOSE  -  DEFINE  INITIAL  CONDITION  FREE  SHAPE  DIGGES  MODEL 
C  -  FOR  ACLS  TRUNKS 

C 

C  HEWTflD  -  SEE  AFFDL-TR-7l^S0.  THEORY  OF  AN  ACLS  FOR  AIRCRAFT 
C  -  ST  KENNEflLT  H.  DIGGES.  JUNE  1371. 

C 

C  LIMITATIONS  -  USED  CH.Y  BY  ROUTINE  QNWT  IN  SUBROUTINE  IC 
C 

C  WRITTEN  BY  -  J.R.KILNER  AND  H.K.UAHI  OCT.  21  1977 
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C _ 
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Figure  37  Subroutine  "IC"  Flowchart  (Continued) 
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SUBROUTINE  103  (X.H.K.I'l.P) 


IW8SE  -  OEfTNEjNJTm^CONDI-riON  LOWED  SHAPES  DIGGE3  MODEL 

mnm  :  ?fiWraKa%"!.“3  r"  WRCWT 

LIMITATIONS  -  USED  OH.T  BT  ROUTINE  QNWT  IN  SUBROUTINE  IC 
WITTEN  IT  -  J.R.  KILMER  AND  HJLHAHI  OCT.  21  1277 
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Figure  37  Subroutine  "IC"  Flowchart  (Concluded) 
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FUNCTION  TBLl(PR,AXO»DPR,NPTS»IS,NA) 


Figure  38  Flowchart  for  Function  "T8L1" 


81 


FUNCTION  TBL2  (PR.ZO.ZORS 


Figure  39  Flowchart  for  Function  "TBL2 


la)  FREE 
TRUNK 
SHAPE 


(b)  LOADED 
TRUNK 
SHAPE 


T 


t 


(c)  TRUNK 

FREE -BODY 
DIAGRAMS 


Membrane  Tr  nk  Element  Geometry  and 
Free- Body  Diagrams 


or 


Tt  ■  <VPc>r2 

similarly  T  =  ^Pt'Pa^rl 
Equating  the  above  two  equations 

rl  Pc  "Pa 

r2  r  r  Pt"Pa 

From  the  membrane  geometry  diagram 

tQ  *  +  r2^2  w>1ere  L0  =  membrane  peripheral  length 

Z0  *  f‘1  (1-cosij)  -  B 
Zo  3  r2  d-cos^) 

A  *  r1  sin^  +  r  sin$2 
Y0  ■  r2  S1n*j 

The  Free  Shape  equation  summary  and  problem  statement  are: 
Solve  iteratively  the  following  four  equations. 


Input  constants: 

Unknowns: 

Determine: 


rl*  r2’*r*2 

Zo  *  r2(l-cos^) 

Y0  -  r2sin^ 

L1  ’  pl*l 

for  0  <  P  <  1.0 
—  r  — 
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The  derivation  for  the  loaded  trunk  shape  is  similar  to  the  free  shape  (see 
Figure  40b).  The  membrane  free-body  is  identical  to  the  free  shape  since  no 
ground  tractions  are  assumed.  Therefore, 


From  the  membrane  geometry  diagram 

l0  ■  *  r^2  ♦  L3 

Z0  ■  r^l-cos^)  +  B 
ZQ  3  r2(l-cos^2) 

L3  +  r2  si  n$2  3  A  -  r^  sin^ 

The  Loaded  Shape  equation  summary  and  problem  statement  are: 
Solve  iteratively  the  following  five  equations. 


Lo  ■  rl*l  *  r2*Z  *  4 
A  3  r1  si  n^  +  r 2  sin$2  +  L3 

B  3  r^l-cos^)  -  r2(l-cos^2) 

Zo  *  ^(i-cos^) 


Input  Constants: 

Unknowns: 

Oetermine: 


rr  r2’  *1*  *2*  L3 
Yq  *  r2  sin*2 
L1  “  rl*l 
L3 

for  0  <  Pr  <  1.0 
and  0  <  ZQ/Zofs  <  1.0 
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where  ZQfs  is  the  free  shape  ZQ  and  can  be  considered  the  point  at  which  the 
membrane  first  makes  contact  with  the  ground  for  a  given  value  of  P^,. 

Two  trunk  section  areas  are  required  for  model  computations.  They  are 


1.  A$  *  trunk  section  area  used  to  calculate  trunk  volume,  and 

2.  Acy  ■  trunk  section  area  inboard  of  ground  contact  point  used 

in  the  cushion  volume  calculation. 


Both  these  areas  can  be  explicitly  determined  from  the  membrane  equation  set 
solution  (i.e.,  r^  r2,  <$2  and  L3) 


As  a  J  r\  (2*i  +  sin^)  +  -[•  r22 

+  (Zo  ’  ri)(A-Yo)  +  riL3 
This  equation  applies  to  both  free 


(2*2-sin 2*2 ) 

and  loaded  trunk  shapes. 
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4.2. 1.2  End  Element 


Two  end  element  models  are  provided  as  a  user  option:  the  constrained 
membrane  and  frozen  model . 

4. 2. 1.2.1  Contrained  Membrane  Model 

The  contrained  membrane  model  is  essentially  the  same  as  the  membrane  side 
element  model  except  no  outward  bulging  occurs  as  the  result  of  trunk  and 
cushion  pressure  variations.  This  is  achieved  using  the  side  element  model 
and  arbitrarily  setting 


4.2. 1.2.2  Frozen  Model 

The  frozen  model  represents  the  trunk  element  as  a  rigid  circular  cylinder 
which  crushes  at  ground  contact  (see  Figure  41). 

Assumptions  1,2, 3, 5, 6  and  9  for  the  side  element  model  also  apply  for  the 
frozen  model.  The  following  set  of  equations  derived  for  the  frozen  model  can 
be  solved  explicitly  and  do  not  require  the  parametric  method  used  for  the 
membrane  model.  However,  for  convenience  within  the  trunk  component  logic, 
the  equations  are  evaluated  for  n  values  of  ZQ  at  T=0  and  the  output 
parameters  stored  in  table  look-up  arrays  similar  to  those  used  for  the 
membrane  model . 


The  equations  for  the  trunk  geometry  as  shown  in  Figure  41  are 


*  r2  sin$2 
C2  3  ~r2  cos<*2 


cos 


Variables  r^,  r 2,  and$2  are  determined  from  the  equations  for  the 

free-shape  membrane  model  for  P^  *  0. 
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■-  'v  ■wmj" 


wxtmjt)  vmsym* 


Then: 


4  •  ri  tW 

L3  -  2rt  Sln01 

4  -20ri 

A  cy-  I  ( V9^  r22  -  T  (Zo’C2)L3  +  7  C2  (CrL3} 

A  s  •  \  CCXC2  ♦  (♦j*Q1)  rx2  +  (^-0j)  r22  +  (Z^tyLj 
♦  (A-C1)(C2-8)a 

The  differences  between  the  constrained  membrane  and  frozen  end  element  models 
are  illustrated  in  Figure  42.  The  frozen  model  predicts  considerably  more 
ground  contact  area  than  the  constrained  model  for  specified  vertical 
deflections.  Increased  contact  area  will  contribute  to: 

1.  Increased  ground  reaction 

2.  Increased  trunk  damping 

3.  Increased  ground-trunk  friction,  and 

4.  Reduced  number  of  trunk  perforations  exposed  to  the  cushion 
and  atmosphere  if  lubrication  is  present. 

4.2.2  Trunk-Cushion  component  Structure 

The  component  equations  and  logic  can  be  best  understood  following  the 
sequence  of  events  within  the  component  computer  code.  As  shown  in  component 
flow  chart  (Figures  35  and  36)  the  calculations  of  model  constants  and  trunk 
element  parameters  are  carried  out  if  T=0  and  the  input  data  differs  from  the 
previous  case.  Area,  volume  and  force  parameters  are  then  set  to  zero  for  the 
suwnation  procedure  in  the  following  loop  in  which  the  element  equations  are 
evaluated  for  each  element  in  the  trunk.  For  each  pass  through  this  loop  the 
following  operations  are  performed: 
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CONSTRAINED 

MEMBRANE 

MODEL 


FROZEN 

MODEL 


4 


TRUNK 

FOOT  PRINT  2 
(INCHES) 

1 
0 

0  12  3 

TRUNK  DEFLECTION 
(INCHES) 


Figure  42  Comparison  of  Constrained  Membrane  and 
Frozen  End  Elements 
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1.  A  test  Is  made  for  side  or  end  element  determination. 


P  »J 
r  1 


P  -P 
c  a 


t  a 


(side  element)  If  B i  *  0 


[  0  (end  element)  if  6  ^  *  0 


2.  The  approximate  location  of  point  “t"  is  determined  in  the  body  axis. 

(See  Figures  33  and  43).  First  the  free-shape  Zo  and  Yo  are  calculated  by 
interpolating  element  parameter  arrays 


Zof$u  ■ 

Zof$  *  Si  Zofsu 

Yofs  ■  S1  f2  (pr‘  W 

where  subscripts  "fs"  denotes  free-shape  and  subscript  "u" 
denotes  unsealed  variables  where  Si  is  a  user  supplied  scaling 
value.  The  functions  f^  refer  the  interpolation  subroutines  for  element 
parameter  data  (see  Figures  38  and  39). 

Then 


*Bt  *  XBa1  *  Yofs  s,n8i 
YBt  ’  a(YSa1  *  Yofs  c°s81> 

ZBt  '  ZBai  *  zofs 

(+1  for  elements  on  right  side 
-1  for  elements  on  left  side 

3.  The  location  of  point  "t"  is  then  determined  in  earth  axis. 


V 

2 

'  X  ’ 

*  [A] 

V 

ZEt 

►  * 

-ALT 

»  * 

YBt 

k  4 

it  4 

.ZBt. 
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Figure  43  Trunk  Geometry 


92 


i 


where  matrix  A  is  the  8ody-to-Earth  Euler  angle  transformation 
(see  Figure  44). 
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and 


[A]  -  [A*]  [Ad]  [A4] 

4.  The  terrain  height  at  location  ”t"  is  evaluated  from  subroutine  TERRA. 

Note  that  for  the  Body  and  Earth  axes  positive  Z  is  down  and  for  the 
terrain  model  positive  Z  is  up  (see  figure  45). 

The  terrain  model  provides  an  option  for  defining  the  surface  profile  with 
a  sinusoidal  function  or  with  a  user  supplied  table  of  profile  elevations. 
These  methods  have  the  following  features: 

a.  Profile  Defined  In  Tabular  Form  (see  Figure  46) 

The  profile  data  must  be  specified  at  even  Increments  along  the 
surface.  Linear  interpolation  Is  performed  between  points.  The 
elevation  outside  of  the  defined  region  takes  the  value  of  the  nearest 
point  (see  flow  chart  in  Figure  47), 
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PROFILE  DEFINED  IN  TABULAR  FORM 


DISCRETE  (1-COS)  OR  CONTINUOUS  SINUSOIDAL  PROFILE 


Figure  46  Terrain  Geometry 
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b.  Discrete  (1-cosine)  or  Continuous  Sinusoidal  Profile 

Any  number  of  sequential  (1 -cosine)  bumps  or  dips  may  be  specified. 

The  profile  begins  at  the  Earth  axis  origin.  A  continuous  sinusoid  is 
represented  by  a  large  number  of  (1-cosine)  bumps. 

Thus 


ZEg  -  TERRA  (XEt) 

5.  A  test  is  made  to  determine  if  trunk  element  is  in  contact  with  the  ground 
(see  Figure  45). 

Zgap  *  '  ZEg  "  ZEt 
Zo  “  Zofs  +  Zgap 

If  Zgap  0  Free  Shape  trunk  calculations  are  performed.  If  not.  Loaded 
Trunk  equations  are  evaluated. 

a.  Loaded  trunk  calculations 


Scaled  trunk  section  properties  are  evaluated. 


Ji  '2  v  ou’  r 
Sf  f3  <Z„U>  Pr> 
S1  U  <Z,u-  Pr> 

S?  f5  'Zou-  Pr> 


If  trunk  element  is  the  frozen  type 

4  "  Si  fl  (Zou*  Pr> 
otherwi se 
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It  should  be  noted  that  for  the  purpose  of  calculating  element  memDrane 
geometry  when  in  ground  contact,  it  is  assumed  that  the  surface  terrain 
under  the  element  is  parallel  to  the  horizontal  trunk  axis.  The 
resulting  error  introduced  is  minor  for  vehicle  roll  angles  less  than 
10°. 


The  trunk  ground  contact  reaction  is  determined  by 


Ft  "  (pfpa>  At 

Next,  the  velocity  of  point  "t"  in  Earth  axis  is 
determined.  The  velocity  components  in  the  ground  plane 
are  rotated  coincident  with  the  body  axis  to  facilitate  the 
ground-trunk  friction  calculation. 


XBt  *  X8ai  *  <V  ?  41  sl"8i 

YUt'  °CX8ai  *  <WL3> 

ZBt  *  ZBai  +  Zo 


Then: 


and 


•  ^ 
A 

.  ' 

*  ^ 

*Et 

( 

u 

P 

YEt 

’  *  [A 8]  [A«  • 

V 

'  +  [B]  - 

q 

ZEt 
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k  4 

>  . 
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The  instantaneous  coefficients  of  trunk-ground  friction  are  then 
determined. 


Uto  =  f  <VEt} 


where  the  function  "f"  is  the  mu-velocity  relationship  for  the  element 


and 


XEt 

A 

YEt 

uty  pto  7^ 

>  ’  -“tx  h 

Ffy  '  ~uty  7t 


The  element  damping  force  is  evaluated  from  the  following  equation. 


Fd  *  °MP  *EtL3Di 

This  relation  for  element  damping  is  based  on  the  following  assumptions. 

1)  Trunk  element  damping  is  due  to  energy  dissipated  by  the  deformation  of 
the  visco-elastic  trunk  material  and  the  trunk-ground  scrubbing  due  to 
vertical  motion  of  the  vehicle.  It  is  assumed  that  this  damping  force  is 
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proportional  to  the  relative  vertical  velocity  between  the  vehicle  hard 
structure  and  the  ground.  Since  ground  rate  of  change  at  point  t  is  not 
available  within  the  model  the  damping  is  assumed  to  be  proportional  to 
the  absolute  vertical  velocity  at  point  t‘  (see  Figure  43). 


2)  It  is  assumed  that  the  damping  force  is  proportional  to  the  footprint 
length  of  the  trunk  element.  This  assumption  can  be  supported  by  the 
argument  that  at  the  instant  of  initial  trunk-ground  contact  the  damping 
force  must  be  zero.  Whereas  at  maximum  trunk  deflection  the  trunk  strain 
rate  and  ground  scrubbing  must  also  be  maximum  for  constant  values  of 

vehicle  vertical  velocity.  An  appropriate  measure  of  these  factors  is  L3. 

3)  It  is  self  evident  that  the  damping  force  must  be  proportional  to  the 
width  of  the  element,  D^. 

The  damping  constant  of  proportionality. 


is  then  a  function  of  the  retarding  force  per  unit  velocity  per  unit  trunk 
contact  area  and  is  therefore  independent  of  trunk  size. 


100 


The  ground  reaction,  friction  and  damping  forces  for  the  element  are 
sunned  with  the  vehicle  forces  and  moments. 

Fx-Fx+Ifx 

Wlfy  . 

pz  •  Fz  -  Ft  -  Fd 

Tx  ■  Tx  -  Zt  +  ^  YBt  -  Ffy  ZBt 

Ty  ’  Ty  +  lFt  *  Fd}  it  +  Ffx  ZBt 

TZ  3  Tz  +  ?fy  XBt  *  Ffx  YBt 

b.  Free  Shape  Trunk  Calculations 

The  trunk-ground  gap  area  is  summed  to  the  total  and  the  scaled  trunk 


secti 

on 

properties 

are 

A 

*  A 

+  Z 

D. 

gap 

gap 

gap 

Yo  3 

Si 

f2 

^Zofsu’ 

Pr> 

L1  a 

Si 

f3 

^ofsu* 

Pr) 

A  ■ 

s 

S? 

f5 

(Zofsu’ 

Pr> 

Acv  ■  s? 

f6 

^  Zofsu 

■  V 

6.  The  element  cushion  area  and  cushion  center  of  pressure  are  determined  and 
the  resulting  element  cushion  reaction  evaluated. 


f°1  <Yai  *  V 


if 0i  =  0 


(Y  R  ,) 


where: 

R 


ai  cos  g  i 


XBc’XBa1  *  <7  Yo  *  7  Rai> 
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The  element  trunk  and  cushion  volumes  are  then  calculated  and  summed  to 
the  totals. 


if  8  i  =  0 
if  8  ^  *  0 


A  test  is  made  to  determine  if  the  trunk  element  has  lubrication  holes. 

If  so,  the  orifice  areas  contributing  to  trunk-to-atmosphere  and 
trunk-to-cushion  flows  are  calculated.  The  flow  chart  for  the  equations 
and  logic  used  in  these  calculations  are  shown  in  Figure  48.  The 
assumptions  and  equations  for  the  trunk  out  flow  model  are  presented  in 
the  following  text. 


1  ■  1th  element 

*  width  of  perforated  area 
■  length  of  trunk  element 

Lpj  ■  Peripheral  distance  from  outboard  attach  point 
(b)  to  edge  of  perforated  area. 

Ap.j  ■  Orifice  area  per  unit  area  of  perforated  surface 
(or  porosity) 

Assumptions: 

a.  Trunk  out-flow  through  perforations  not  in  contact  with  ground  is  a 
function  of  Pt,  Pc,  Pfl,  Cg1  and  the  corresponding  portion  of  orifice 
area. 

b.  Outflow  from  the  flattened  portion  of  trunk  is  split  based  on  the 
location  of  a  stagnation  pressure  and  an  assumed  triangular 
pressure  gradient. 


c.  The  stagnation  pressure  in  the  gap  is  assumed  to  be  equal  to  the  trunk 
pressure. 

d.  The  gap  pressure  gradients  on  either  side  of  the  stagnation  point  are 
assumed  to  be  equal  (see  Figure  49),  i.e. , 
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Therefore: 


VPc  L3c 
Pfpa  *  L3a 


L3a  +  L3c 


or 


(Pt  -  pa> 

L3a  ’  Oyp.Mfyr)-  L3 


<prpc> 


l3c  (vpa)-(pt-pcr  4 


Consider  the  control  volume  shown  in  Figure  49. 
dQ  K 


Q+dQ 


STAGNATION  POINT 


dl 


dQ  .  0,  Ap1  C02  / 1  (Pt-P(l)) 


Pt*p, 

where  P ( 1 )  s  Pt  — i -  1  for  the  atmosphere  side. 

1  L3a 

Therefore,  the  total  flow  to  the  atmosphere  from  the 
flattened  area  is 


1.2  ■/  "'“I  Api  CD2  4  (VPa>t^  d' 


or 


<>a2  '  !  4a  Ap1  4  c02$  <VPa> 
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Pneumatic  Outflow  and  Pressure  Gradients  for  Flattened  Trunk 


and  similarly 


°C2  ■  T  4c  Ap1  ®1  'WS  (%  -  Pc> 

Thus,  the  flow  from  the  flattened  area  Lja  Di  (or  Ljc  D^)  is  simply  the 
normal  discharge  affected  by  Pt-Pa  (or  Pt-pc)  reduced  by  a  factor  of  2/3. 

Now,  consider  the  condition  in  which  the  perforated  area  makes  up  a 
fraction  of  the  flattened  area.  The  diagram  in  Figure  50  shows  the 
possible  arrangements  of  perforated  area  to  flattened  area. 

Conditions  1,4,10,11,12,  and  13  are  modeled  precisely.  The  area 
relationships  for  conditions  2,3, 5,6, 7,8  and  9  are  accurately  represented 
but  no  adjustment  has  been  made  to  the  triangular  pressure  gradient  for 
flow  in  the  flattened  portion  of  the  trunk. 


It  should  be  noted  that 
section. 


Q_  ,  U  ,  and  L-,  are  zero  for  a  free  truck 
c2  Ja  Jc 


Flow  Area  Adjustments  for  End  Elements. 


Ea2  *  (RN  +  L3  "  7  L3a)  /R0 
Ec2  *  (RN  +  7  L3c)  /R0 
Ecl  *  (RN  '  7  LhclJ  /RD 
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LOADED  TRUNK 


h»«i  ’  Lr  LP 

Lhcl  ■  Lp*  **h”  Ll“  L3 


ATM 

OUTB* D  SIDE 
ATTACH 


FLATTENED  CUSHION 

AREA  SIDE 


SEE  IMS 


Lhal  *  Ll'  Lp 
Lhcl  ■  Lp+  V  L1 


ATM. 

SIDE 


Figure  50  Schematic  of  Possible  Arrangements  of  Trunk 
Perforated  Area  to  Flattened  Area  (Concluded) 
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where  Eflj,  E^,  E^  and  Ecj  are  the  area  adjustment  factors  for  areas  al, 
a2,  c2  and  cl  respectively. 

And  R  »  R  .  +  Y 
N  ai  o 

R0  *  Rai  +  Yofs 

8.  After  perforation  area  is  calculated  the  model  indexes  to  the  next  element 
and  repeats  all  calculations  from  paragraph  1,  Section  4.2.2  in  this 
test.  Once  all  elements  have  been  evaluated  the  trunk  and  cushion  volume 
rates  are  calculated  through  integration  using  a  first  order  lag. 

Vc  -i  <v;-vc) 

vt  -i  (V*  -  Vt) 
where: 

V*,  V*  »  calculated  volumes 

Vc,  Vt  *  volume  state  variables 

Vc,  Vt  *  volume  rate  of  change 
T  »  time  constant 

9.  The  total  effective  areas  for  air  flow  from  the  trunk  to  the  cushion, 
trunk  to  the  atmosphere  and  cushion  to  the  atmosphere  are 
calculated. 

^relief  '  f  <%  ' 

^ca  "  CDgap  Agap 

^ta  "  CD1  Ahal  *  I  C02  Aha2  +  CDa  belief 

Wtc  "  C01  Ahcl  *1  C02  Ahc2 
where  »  trunk  to  atmosphere  relief  valve  area. 


Ill 


The  relief  valve  area  Is  a  function  of  the  pressure  difference  pt_pa  and 
Is  specified  by  the  user. 

10.  Mass  flow  rates  are  calculated  from  subroutine  FNFLOW  using  the  Chester 
Smith  compressible  flow  function  given  the  upstream  and  downstream 
pressure,  temperature  and  effective  area. 


60  CA 


ca 


60  CA. 


tc 


-  60  CAU 


11.  Trunk  and  cushion  pressure  rates  are  evaluated  from  the  following 

functions  derived  from  the  conservation  of  energy  of  the  trunk  and  cushion 
vol umes , 


t 


c 


Where  k  Is  the  polytropic  constant  for  air. 


112 


SECTION  V 


BOEING  ELASTIC  TRUNK  AND  CUSHION  COMPONENT 


5.1  General  Model  Description 

5.1.1  Trunk  Geometry 

The  model  represents  the  geometry  and  pneumatic  characteristics  of  an  elastic 
trunk  enclosing  a  cushion  of  air  (Figure  29). 

The  trunk  is  described  by  a  finite  number  of  elements  which  need  not  be 
uniform  in  cross  section,  width,  or  location  on  the  vehicle  (Figure  30). 

The  trunk  is  composed  of  two  basic  element  types:  side  elements  and  end 
elements  (both  fore  and  aft).  The  side  elements  have  unconstrained  outward 
movement  subject  to  ground  contact  or  changes  in  the  cushion/trunk  pressure 
ratio  and  are  represented  with  a  two-dimensional  membrane  model  as  shown  in 
Figure  31.  The  end  elements  experience  little  outward  motion  due  to 
constraining  peripheral  (hoop)  stresses  and  are  represented  by  a 
three-dimensional  membrane  model  as  shown  in  Figure  51. 

The  element  locations  are  specified  in  the  trunk  coordinate  system  by  the 
position  of  the  trunk  inboard  attach  point  (Figure  33).  The  trunk  axis  is 
referenced  to  the  vehicle  body  axis  as  shown  in  Figure  29.  Thus  the  entire 
trunk  assembly  can  be  repositioned  on  the  vehicle  by  changing  the  value  of  two 
parameters,  BST  and  WLT. 

5.1.2  Pneumatic  System 

A  schematic  of  the  trunk-cushion  pneumatic  system  is  shown  in  Figure  34. 
Compressible  flow  equations  are  used  to  describe  the  state  of  the  fluid.  The 


Figure  51  End  Element  Membrane  Model 


remaining  air  supply  system  not  shown,  such  as  bleed  air,  ejectors,  fans,  etc. 
are  modeled  with  other  EASY  components.  Air  can  be  supplied  to  the  trunk 
and/or  cushion.  Air  outflow  occurs  through  the  trunk  perforations ,  the 
trunk-ground  gap  and  the  trunk  relief  valve. 

5.1.3  Component  Structure 

A  simplified  flow  chart  of  the  trunk  component  is  shown  in  Figure  35  which 
describes  the  basic  function  of  the  model.  Input  variables  required  are 
velocities  and  displacements  from  the  vehicle  rigid  body  component  and 
flow-rate  and  temperature  from  the  air  supply  components.  During  the  initial 
pass  through  the  component  (i.e.  when  T=0)  parameter  data  tables  are 
calculated  from  a  set  of  implicit  equations  which  represent  the  trunk  section 
properties.  Subsequent  evaluation  of  the  element  section  properties  can  then 
be  made  economically  by  interpolating  these  tabular  data.  A  "DO"  loop  is 
entered  within  the  code  which  evaluates  the  individual  element  properties  and 
calculates  their  accumulative  effect  as  a  trunk-cushion  system.  Note  that  the 

terrain  model  is  called  from  the  trunk-cushion  model  and  is  not  an  independent 
EASY  component. 

Vehicle  forces  and  moments  are  calculated  as  well  as  pressure  and  volume  rates 
of  change  and  flow  rates  which  are  outputted  to  their  corresponding  components. 

5.2  Discussion  of  Model  Idealization 

A  detailed  flowchart  of  the  trunk-cushion  component  is  shown  in  Volume  II, 

Table  115  corresponds  to  the  simplified  chart  of  the  model  shown  in  Figure  35. 
The  detailed  description  of  the  approach,  assumptions  and  equations  of  the 
model  subcomponents  will  be  presented  in  the  same  order  as  they  are  shown  on 
the  flowchart. 
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5.2.1  Trunk  Section  Properties 


Each  element  of  the  trunk  is  represented  by  an  elastic  membrane  in  the.shape 
of  either  a  cylinder  (for  side  elements)  or  a  torus  section  (for  end 
elements).  The  sets  of  equations  which  describe  these  shapes  are  implicit  and 
require  an  iterative  solution.  Solving  these  equations  by  iterating  for  each 
call  from  the  EASY  integrator  would  be  very  time  consuming  and  impractical. 

The  approach  used  for  this  model  is  to  initially  solve  the  sets  of  equations 
(i.e.  at  TIME=0)  and  store  the  parameter  data  in  arrays.  The  section 
properties  are  then  evaluated  from  the  arrays  using  table  look-up  routines 
during  model  execution. 

The  section  parameter  data  is  generated  by  subroutine  ELAS  (Volume  II,  Section 
3.7.12).  The  model  will  generate  and  store  up  to  eight  sets  of  parameter  data 
each  representing  an  element  of  size  and  shape  specified  by  the  user. 


5.2. 1.1  Side  Element  (Two-Dimensional  Membrane  Model) 

The  free-body  diagrams  of  a  trunk  side  element  membrane  and  the  trunk  geometry 
for  loaded  and  free  shape  are  shown  in  Figure  40.  The  following  assumptions 
were  observed  in  the  derivation  of  the  governing  equations. 

1.  The  trunk  element  behaves  as  a  two-dimensional  elastic  membrane. 

2.  Trunk  membrane  tension  is  constant  along  its  minor 
periphery  (i.e.,  no  contribution  from  ground  tractions, 

5«.c . ) 

3.  Pressure  change  from  P  to  P  occurs  at  the  venturi 

C  3 

throat. 

4.  The  trunk  membrane  is  tangent  at  ventumi  throat. 

5.  Perforation  nozzle  reactions  are  negligible. 

6.  The  pressures  on  both  sides  of  membrane  element  in 
contact  with  the  ground  are  equal  to  P^. 
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7.  The  portion  of  membrane  in  contact  with  the  ground  is  a 
straight  line. 

8.  The  trunk  is  tangent  at  the  edge  of  ground  contact. 

9.  The  pressure  change  from  Pc  to  P„  and  from  P^  to  Pa 
occurs  instantaneously  at  the  edge  of  the  flattened 

area. 


The  following  derivation  is  for  the  free  trunk  shape  shown  in  Figure  40a. 


From  the  membrane  free  body  diagram 


2Ttsin  2^Z~  2  ^Pt"PJ  r2  sin  2  ^  2 


or 


Tt  3  (,VPc)r2 

similarly  T  =  ^Pt"Pa^rl 


Equating  the  above  two  equations 

rl  pc  -Pa 

J.  .  1  -  ?  where  P  -  fa-i 

r2  r  t  a 


From  a  consideration  of  the  membra-ne  elastic  characteristics,  the  meridian 
strain  is  defined  by  a  user-input  table  to  be  Z<p  =  f(Tt).  The  membrane 
meridian  length  L  may  be  represented  as:  % 

» 

L  -  Lf (i  +  e^)  ’ 

where 

W  ■  V'1  *  e*i>  • 


1^  ■  mentrane  meridian  free  length 
Lq  •  membrane  meridian  installed  length 

■  meridian  strain  in  the  installed  configuration  (unpressuri zed) . 
C+  9  meridian  strain. 


117 


From  the  membrane  geometry  diagram 


L  *  ri*  l  +  r2*  2 
Zg  *  Tj  (l-cos^)  -  8 

Zo  3  r2  ^•cos^2) 

A  =  r1  s i +  r?  sin<$2 

^o  *  r2  s^n^2 

The  Free  Shape  equation  summary  and  problem  statement  are: 
Solve  iteratively  the  following  five  equations. 


L  =  rl*l  +  r2^2 
A  «  Tj  sin<^  +  r2  sin$  2 

B  =  fj  (l-cos<frj)  -  r^  (l-cos<$2 ) 

Input  constants:  Ppf  A,  B,  LQ,  £<pi 

Unknowns:  rl •  r2*  *1 *  $2'  e<P 

Determine:  ZQ  «  ^(l-cosL) 

Y0  »  r2sin^2 

4  *  ri*i 

for  0  _<  ?r  <_  1.0 

The  derivation  for  the  loaded  trunk  shape  is  similar  to  the  free  shape  (see 
Figure  40b).  The  membrane  free-body  is  identical  to  the  free  shape  since  no 
ground  tractions  are  assumed.  Therefore, 
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From  the  membrane  elastic  characteristics 


U 

L 


From  the  membrane  geometry  diagram 


ZQ  *  r^(l-cos^)  +  B 

» 

ZQ  =  r2(l-cos*2) 

+  r2  sin$„  *  A  -  r  sin^ 

The  Loaded  Shape  equation  sunmary  and  problem  statement  are: 
Solve  iteratively  the  following  six  equations. 


L 

L 


A  *  r^  sin^1  +  r 2  sin$2  + 

B  *  r^l-cos^)  -  r2(l-cos«2) 
Z0  *  r2(l-cos*2) 
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Input  Constants: 

Unknowns: 

(Determine: 


pr»  A*  0.  z0»  LQ,  Eei 
rl*  r2’  *1’  *2*  L3’ 

\  ■  r2  s,n4 

4  ’  44 

L3 

for  0  £  P  £1.0 

o"<  Zo/ZQfs  <  1.0 


where  ZQfs  is  the  free  shape  ZQ  and  can  be  considered  the  point  at  which  the 
membrane  first  makes  contact  with  the  ground  for  a  given  value  of  P  . 

J  p 

Two  trunk  section  areas  are  required  for  model  computations.  They  are 


1*  As  *  trunk  section  area  used  to  calculate  trunk  volume,  and 
2*  Acv  3  trunk  section  area  inboard  of  ground  contact  point  used 
in  the  cushion  volume  calculation. 


Both  these  areas  can  be  explicitly  determined  from  the  membrane  equation  set 
solution  (i.e.,  r^  r2,  ♦j,  *2  and  l3) 

Area  A$ 


As  *  T  r\  (2*1  *  sin2*x  )  +  |  r22  (2*2-sin2*2) 


+  <Zo  ‘  rlKA-Y0)  ♦  rfo 
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This  equation  applies  to  both  free  and  loaded  trunk  shapes. 


Area  A 


cv 


Acv  *  T  ^V*1"2^ 

5.2. 1.2  End  Element 

The  end  element  trunk  geometries  for  loaded  and  free  shape  are  shown  in  Figure 
51.  The  following  assumptions  were  observed  in  the  derivation  of  the 
governing  equations. 

1.  The  trunk  element  behaves  as  a  three-dimensional  elastic  membrane. 

2.  Hooke's  law  is  applicable  for  the  consideration  of  combined  loads 
and  strains. 

3.  A  membrane  meridian  is  comprised  of  two  connected  elliptical  arcs. 

4.  Pressure  change  from  Pc  to  occurs  at  the  venturi  throat. 

5.  A  trunk  meridian  is  tangent  at  the  venturi  throat. 

6.  Perforation  nozzle  reactions  are  negligible. 

7.  Pressures  on  both  sides  of  the  membrane  element  in  contact  with  the 
ground  are  equal  to  Pt« 

8.  The  portion  of  the  membrane  meridian  in  contact  with  the  ground  is  a 
straight  line. 

9.  A  trunk  meridian  is  tangent  at  the  edge  of  ground  contact. 

10.  The  pressure  change  from  Pc  to  and  from  Pt  to  Pa  occurs 

instantaneously  at  the  edge  of  the  flattened  area. 
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A  trunk  end  is  assumed  to  be  a  surface  of  revolution  obtained  by  rotation  of  a 
plane  curve  about  an  axis  lying  in  the  plane  of  the  curve.  The  curve  is 
called  the  meridian,  and  its  plane  is  a  meridian  plane.  A  finite  element  of 

the  trunk  surface  is  cut  out  by  two  adjacent  meridians  and  two  parallel 

circles  as  shown  in  Figure  52a. 

The  position  of  a  meridian  is  defined  by  the  angle  0,  measured  from  some  datum 
meridian  plane;  and  the  position  of  a  parallel  circle  is  defined  by  the  angle 
made  by  the  normal  to  the  surface  and  the  axis  of  rotation.  The  meridian 
plane  and  the  plane  perpencicular  to  the  meridian  are  the  planes  of  principal 
curvature,  and  the  corresponding  radii  of  curvature  are  denoted  by  r^  and  r ^ 
respectively.  The  radius  of  a  parallel  circle  is  denoted  by  rQ. 

The  magnitudes  of  the  normal  forces  per  unit  length  are  denoted  by  Ng  (hoop 

load)  and  (meridian  load)  as  shown. 

Equilibrum  Equations  (See  Reference  5) 

An  equation  of  equilibrium  is  obtained  by  summing  up  the  projections  of  the 
forces  in  the  direction  normal  to  the  surface  of  the  element  shown  in  Figure 
52a.  The  forces  acting  on  the  upper  and  lower  sides  of  the  element  have  a 
resultant  in  the  normal  direction  equal  to  N$-ro.d0.d$  . 

The  forces  acting  on  the  lateral  sides  of  the  element  have  a  resultant  N0-r^-d^-d0 
in  the  radial  direction  of  the  parallel  circle.  The  normal  component  is  equal 
to  Ng.r^«sin$.d$.d0. 

The  pressure  load  acting  normal  to  the  element  is  -P* r^* rQ*d^*d0. 

Sunmation  of  these  forces  yields  an  equation  of  equilibrium 

N*.r0  +  Ng.rj.sin^  -  P-rj*^  *  0. 
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Dividing  by  ry  rQ  and  noting  that  rQ  *  sin*,  yields 


A  second  equation  may  be  derived  by  considering  the  equilibrium  of  the  trunk 
surface  formed  by  the  rotation  of  the  trunk  meridian  about  a  vertical  axis 
(Figure  52b).  Considering  the  equilibrium  of  the  trunk  portion  represented  by 
rotation  of  the  line  ST  yields: 

21T  Rq  sin*  =  IT P(RQ2  -  b2)  or 

PP2 

N*  *  ^  (Ro  +  b)’ 

The  cross  sectional  shape  of  a  loaded  end  element  meridian  is  shown  in  Figure 
53.  The  meridian  shape  is  assumed  to  be  two  elliptical  arcs  connected  by  a 
straight  line  segment  (L3)  which  is  in  contact  with  the  ground. 


r2  P2  + 

rQ  -  r2  sin* 


A  similar  analysis  is  performed  for  the  inner  ellipse  to  determine  equilibrium 
loads  over  the  inner  meridian. 


The  trunk  membrane  is  an  elastic  material  and  as  such  will  experience  a 
deformation  upon  the  application  of  a  load.  This  section  presents  the 
relationships  between  load  and  strain  for  the  trunk  membrane  and  the 
associated  assumptions. 
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Figure  53 


AXIS  OF 
REVOLUTION 


End  Element  Geometry 


'  / 
V 


Figure  54  uninflated  End  Element  with  Radial  Construction 


Strain  €  ,  is  defined  as 


L  -  L  L 

€  ■  — j -  ■  7—  -1,  where 

L  ■  Instantaneous  length 
Lf  ■  free  length  (no  loading). 

Combined  loads  and  strains  in  the  meridian  (*)  and  hoop  (8)  directions  are 
handled  as  follows: 


1.  Oirectional  loads  and  Ng  are  calculated  from  the  equilibrium 
equations; 

2.  Combined  loadings  are  then  defined  by 


V  -  N*  -  VNg 


N0  -  Vl% 


where 


V  is  Poissons  ratio,  a  user  Input. 

3.  The  combined  loads  N$c  and  Nqc  are  used  to  determine  the  meridian 
and  hoop  strains  from  the  user- input  load  deflection  curves 

£»-  U  (N*c) 
ee"  fe  (N0c ) 

If  the  membrane  for  the  trunk  ends  is  constructed  of  radial  cords  which  bear 
most  of  the  meridian  load  (rather  than  the  material  between  cords),  a 
modification  to  the  previous  load/deflection  equations  must  be  made. 

Figure  54  depicts  a  section  of  end  membrane  with  radial  cords  with  a  meridian 
force  T  applied.  The  meridian  loads  at  the  outer  (r2)  and  inner  (r^)  radii 
are  respectively 

■  T/^0  (5-1) 


N*2  -  T/r2Q 


(5-2) 
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Hooke's  law  gives: 


e,  ■  jj-  in,  -  vn9)  .  ,  die  ■  o)  (5'3) 

where  is  a  material  constant. 

Substitution  of  equations  (5-1)  and  (5-2)  into  (5-3)  would  yield  differing 
strains  at  the  inner  and  outer  radii.  This  is  contrary  to  the  assumption  that 
the  radial  cords  assumed  the  load,  resulting  In  a  constant  tension  force  and 
therefore,  a  constant  strain  along  the  meridian  length  of  the  cords. 

Thus,  for  radial  cord  construction,  equation  (5-3)  is  modified  to  yield 


where 


(5-4) 


r$  is  a  reference  radius  such  that  at  r  ■  r$,  the  discrete  membrane 
cord  construction  can  be  treated  as  a  continuous  homogeneous 
material . 


Equations  of  Geometric  Compatibility 

Previous  sections  have  discussed  the  calculation  of  loads  and  subsequent 
strains  for  the  elastic  membrane.  The  final  step  In  the  analysis  of  the  trunk 
membrane  is  a  consideration  of  the  membrane  deformation  and  its  associated 
conditions  of  geometric  fit.  That  is,  the  deformations  which  result  from  the 
calculated  loads  must  yield  a  membrane  shape  which  In  turn  will  yield  the  same 
loads. 

In  the  meridian  direction  the  membrane  length  can  be  divided  into  three 
distinct  segments  (Figure  55). 


W 


^3T 


A  -  OUTER  ATTACH  POINT 
B  -  INNER  ATTACH  POINT 


UNINFLATEQ 

MEMBRANE 


EDGE  VIEW  OF 
UNINFLATED  MEMBRANE 


Figurt  55  End  Eltment  Geomttrlc  Compatibility 


Lj  peripheral  distance  from  outer  attach  point  to  ground  contact; 

L2  peripheral  distance  from  inner  attach  point  to  ground  contact; 

ground  contact  length. 


For  the  inflated  membrane: 


TT 


e*o 


-1 

+ 


u2 

TT 


*2 


1  +  £^2 


*  meridian  free  length 
Lq  =  meridian  installed  length  (uninflated) 

£^o  =  initial  meridian  strain  in  uninflated  configuration 

£$^  =  calculated  meridian  strain  over 

£+2  3  calculated  meridian  strain  over  l_2 

£$3  =  calculated  meridian  strain  over 


The  meridian  lengths  and  are  elliptical  arcs  (by  assumption)  which  can 
be  determined  by  evaluation  of  an  elliptic  integral  of  the  second  kind. 


The  strains  £^  and  £ ^  are  found  by  calculating  the  average  meridian  loads 
over  and  using  the  equilibrum  equations  and  the  user-defined 
load/deflection  tables.  For  the  loaded  shape  £^  is  defined  as  the  average 

£  *3  3  *5  (£$j  +£(>2) 


In  the  hoop  direction,  the  strain  is  defined  by  considering  a  membrane  section 
of: 


ce 


rA0-  rfA0 
TfAQ 


where 


(5-5) 
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r^  *  radius  to  a  Doint  on  the  membrane  in  the  free  condition 

r  =  radius  to  the  same  point  in  the  stretched  condition. 

The  hoop  strain  present  in  the  installed  uniflated  membrane  configuration 
is  defined  by: 


’0o  = 


r  -  r, 
o  f 


where 


(5-6) 


rQ  is  the  radius  to  a  point  in  the  installed  uninflated  configuration. 


Substituting  equation  (5-6)  for  r^  into  equation  (5-5)  yields: 


ce  •  -r  (i +  ee0>  -  (  } 

0 

Equation  (5-7)  must  be  satisfied  for  the  membrane  over  the  meridian  lengths 
and  l_2  individually  that  is: 

Eei  'T^7  (1  ‘ 1 

£92  ■  TFTT  »  **%>  -  1 

0  Z 

The  radius  r^  is  determined  by  averaging  the  radii  to  points  on  the  membrane 
surface  over  L^.  Likewise  is  determined  by  averaging  the  radii  to  points 
on  the  membrane  surface  along  L^. 

The  installed  radii  (r  ),  and  (r  )_  are  calculated  from  a  consideration  of  the 

o  1  o  2 

membrane  in  the  installed  uninflated  configuration.  An  edge  view  of  a  trunk 
meridian  length  Lq  may  be  subdivided  into  lengths  L^q  and  L^q  which  are 
defined  as  the  installed  meridian  lengths  which  correspond  to  lengths  ,  I2, 
and  L^. 
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Since: 


4  -  4f 


Lio  "  4f 


it  follows  that: 


4  '  {1  +  £*o} 


l2*(i  +e*0) 

tttt; 


l3  *  (1  +  UJ 

L30  =  1  +£*3 

The  installed  radii  (r  ),  and  (r  )„  are  then  defined  as  the  respective  radii 

Ol  0  2 

to  the  midpoints  of  L^g  and  l_2g. 

Equation  Summary 

The  free  shape  equation  suirmary  (refer  to  Figure  5-3)  and  problem  statement 


Solve  iteratively  the  following  equations: 


4o +  4o 

Yrt  +  za 
o  a 


-  R  cos*  .  -  F  +  R^  cos*  , 
a  1  e  c. 


t^T  (I  *  £«») ' 1 


'v 


E02  '  (l  *  C9o>  '  1 

F  »  PR  A2  V/Q2 


Input  constants:  LQ,  A,  3,  Cqq,  £*0,  Pt,  Pc 
Unknowns:  Q,  V,  ^  E,  F,  *2 

The  loaded  shape  equation  summary  and  problem  statement  are 
Solved  iteratively  the  following  equations: 


Lo  *  L10  +  L20  +  L30 

A  »  Y  +  Z  +  L, 

o  a  o 

B  =>  V  -  R.  cos  $  .  -  F  +  R  cos  $  0 
a  1  e  i. 

£9i  -  ^  »  ♦  £90)  -  1 


Input  constants:  Lq,  A,  B,  £Qo,  e*Q,  ?v  ?q 
Unknowns:  Q,  V,  ^  E,  F,  *2,  L-j 
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5.2.2  Trunk-Cushion  Component  Structure 


The  component  equations  and  logic  can  be  best  understood  following  the 
sequence  of  events  within  the  component  computer  code.  As  shown  in  component 
flow  chart  (Figure  35)  the  calculations  of  model  constants  and  trunk  element 
parameters  are  carried  out  if  T=0  and  the  input  data  differs  from  the  previous 
case.  Area,  volume  and  force  parameters  are  then  set  to  zero  for  the 
summation  procedure  in  the  following  loop  in  which  the  element  equations  are 
evaluated  for  each  element  in  the  trunk.  For  each  pass  through  this  loop  the 
following  operations  are  performed: 

1.  A  test  is  made  for  side  or  end  element  determination. 

2.  The  approximate  location  of  point  "t"  is  determined  in  the  body  axis. 

(See  Figures  33  and  43).  First  the  free- shape  Zo  and  Yo  are  calculated  by 
interpolating  element  parameter  arrays 

Vs  ■  fl  <Pr-  Pt> 

Vs  *  f2  <Pr-  Pt> 


where  subscripts  "fs"  denotes  free-shape.  The  functions  f^  refer  the 
interpolation  subroutines  for  element  parameter  data  (see  Figures  38  and 
39) 

Then: 


XBt  ■  XBai  +  Vs  s1n8i 

V°(Vi  +  Vsc°sl3i) 

ZBt  "  ZBaf  *  Zofs 


where : 


1+1  for  elements  on  right  side 
-1  for  elements  on  left  side 
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3.  The  location  of  point  "t"  is  then  determined  in  earth  axis. 


*  ' 
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where  matrix  A  is  the  Body-to-Earth  Euler  angle  transformation 
(see  Figure  44). 
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4.  The  terrain  height  at  location  "t"  is  evaluated  from  subroutine  TERRA. 
Note  that  for  the  Body  and  Earth  axes  positive  Z  is  down  and  for  the 
terrain  model  positive  Z  is  up  (see  Figure  45). 

The  terrain  model  provides  an  option  for  defining  the  surface  profile 
with  a  sinusoidal  function  or  with  a  user  supplied  table  of  profile 
elevations.  These  methods  have  the  following  features: 
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a.  Profile  Defined  in  Tabular  Form  (see  Figure  46) 

The  profile  data  must  be  specified  at  even  increments  along  the 
surface.  Linear  interpolation  is  performed  between  points.  The 
elevation  outside  of  the  defined  region  takes  the  value  of  the 
nearest  point  (see  flow  chart  in  Figure  47). 

b.  Discrete  (1-cosine)  or  Continuous  Sinusoidal  Profile 

Any  number  of  sequential  (1-cosine)  bumps  or  dips  may  be  specified. 
The  profile  begins  at  the  Earth  axis  origin.  A  continuous  sinusoid 
is  represented  by  a  large  number  of  (1-cosine)  bumps. 

Thus 

ZEg  =  TERRA  (XEt) 

5.  A  test  is  made  to  determine  if  trunk  element  is  in  contact  with  the 
ground  (see  Figure  45). 

Zgap  =  '  ZEg  '  ZEt 

1*1  +  Z 

Lo  ofs  gap 

If  Z  >0  Free  Shape  trunk  calculations  are  performed.  If  not.  Loaded 
gap  — 

Trunk  equations  are  evaluated, 
a.  Loaded  trunk  calculations 

Trunk  section  properties  are  evaluated. 
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It  should  be  noted  that  for  the  purpose  of  calculating  element 
membrane  geometry  when  in  ground  contact,  it  is  assumed  that  the 
surface  terrain  under  the  element  is  parallel  to  the  horizontal 
trunk  axis.  The  resulting  error  introduced  is  minor  for  vehicle 
roll  angles  less  than  10°. 

The  trunk  ground  contact  reaction  is  determined 


-  DiL 
*  A 


Next,  the  velocity  of  point  "t"  in  Earth  axis  is  determined.  The 
velocity  components  in  the  ground  plane  are  rotated  coincident  with 
the  body  axis  to  facilitate  the  ground-trunk  friction  calculation. 


sinfl. 


YBt  *  °  [YBai  +  <Y„  *  7  4>  c05M 
*8t  '  Z8af  *  Zo 


Then: 
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The  instantaneous  coefficients  of  trunk-ground  friction  are  then 
determined. 


f  (V  ) 
Et 


where  the  function  "f"  Is  the  mu-velocity  relationship  for  the 
element  and  is  defined  as  shown  below. 
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The  element  damping  force  is  evaluated  from  the  following  equation. 

Fd  "  °MP  ^EtL3°1 

This  relation  for  element  damping  is  based  on  the  following 

assumptions. 

1)  Trunk  element  damping  is  due  to  energy  dissipated  by  the 
deformation  of  the  visco-elastic  trunk  material  and  the 
trunk-ground  scrubbing  due  to  vertical  motion  of  the  vehicle. 

It  Is  assumed  that  this  damping  force  is  proportional  to  the 
relative  vertical  velocity  between  the  vehicle  hard  structure 
and  the  ground.  Since  ground  rate  of  change  at  point  t  is  not 
available  within  the  model  the  damping  is  assumed  to  be 
proportional  to  the  absolute  vertical  velocity  at  point  t'  (see 
Figure  43). 

2)  It  is  assumed  that  the  damping  force  is  proportional  to  the 
vertical  deflection  of  the  trunk  element.  This  assumption  can 
be  supported  by  the  argument  that  at  the  instant  of  initial 
trunk-ground  contact  the  damping  force  must  be  zero.  Whereas 
at  maximum  trunk  deflection  the  trunk  strain  rate  and  ground 
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scrubbing  must  also  be  maximum  for  constant  values  of  vehicle 
vertical  velocity.  An  appropriate  measure  of  these  factors  is 

L3‘ 

3)  It  is  self  evident  that  the  damping  force  must  be  proportional 
to  the  width  of  the  element,  D^. 

The  damping  constant  of  proportionality. 


is  then  a  function  of  the  retarding  force  per  unit  velocity  per 
unit  trunk  contact  area  and  is  therefore  independent  of  trunk 
size. 


The  ground  reaction,  friction  and  damping  forces  for  the 
element  are  sunned  with  the  vehicle  forces  and  moments. 
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.  Free  Shape  Trunk  Calculations 


Trunk  section  properties  are  evaluated  and  the  trunk-ground  gap  area 
is  sunned  to  the  total. 


A  =  A  +  Z  D . 
gap  gap  gao  i 

Yo  '  f2  <Pr*  Pt> 

L 


'1 


*  f3  (Pr*  V 


13c 


U  (Pr-  Pt> 


Ac,  *  f6  <Pr-  Pt> 


The  element  cushion  area  and  cushion  center  of  pressure  are  determined 
and  the  resulting  element  cushion  reaction  evaluated. 


»1  <V,1  *  V 


1.  <Yo.V 

?  '  (V'of S**a,> 


if  B  4-0 


if B  i  t  0 


where: 


ai  cosg. 


xBc  -  xBai  *  <4  Yo  -  T  Rai>  sinS  , 

Y8c  I  <vBaf  *  V 

YBai  *  <!  Yo  -  7  *a.>  “*1 


if  »  0 
if  ei  *  0 


Fc  ■  <pc-pa>A= 

F2  *  F2  "  FC 

Tx‘Tx-!c  YBc 
Ty-“  Ty  +  Fc  xBc 

The  element  trunk  and  cushion  volumes  are  then  calculated  and  sunned  to 
determine  the  total  volumes. 


V*  -  V*  +  Z  A  -  0.  A  « 
C  C  0  c  1  cv 


2/3T  ♦  R  . 

0  31 

Y  +  R 
'o  ai 


iffl  .  *  0 


v! 


v*t  + 


D,  A 


7.  A  test  is  made  to  determine  if  the  trunk  element  has  lubrication  holes. 

If  so,  the  orifice  areas  contributing  to  trunk-to-atmosphere  and 
trunk-to-cushion  flows  are  calculated.  The  flow  chart  for  the  equations 
and  logic  used  in  these  calculations  are  shown  in  Figure  48.  The 
assumptions  and  equations  for  the  trunk  out  flow  model  are  presented  in 
the  following  text.  Definition  of  trunk  perforations: 


*"hi  *  width  of  perforated  area 
■  length  of  trunk  element 

Lpi  *  Peripheral  distance  from  outboard  attach  point 
(b)  to  edge  of  perforated  area. 

Apf  *  Orifice  area  per  unit  area  of  perforated  surface 
(or  porosity) 

Assumptions: 

a.  Trunk  out-flow  through  perforations  not  in  contact  with  ground 
is  a  function  of  P  P  p  (V  and  the  corresponding  portion 
of  orifice  area. 
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b.  Outflow  from  the  flattened  portion  of  trunk  Is  split  based  on 
the  location  of  a  stagnation  pressure  and  an  assumed  triangular 
pressure  gradient. 

c.  The  stagnation  pressure  In  the  gap  is  assumed  to  be  equal  to  the 
trunk  pressure. 

d.  The  gap  pressure  gradients  on  either  side  of  the  stagnation 
point  are  assumed  to  be  equal  (see  Figure  49),  i.e., 


a  ■  a 


Therefore: 


.  <Vc> 

l3c  (ft-fr,)Tpt'-rc7  4 

Consider  the  control  volume  shown  below; 
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«  ■  °1  »pi  «D2  Vf  (pt  -  P<i>> 


where 


P,-P_ 

P(l)  ■  P^  — ^ — *  1  for  the  atmosphere  side. 
3d 


Therefore,  the  total  flow  to  the  atmosphere  from  the 
flattened  area  is: 

L3a 

%z-fa  °i  V  co2  Vs  (W-qa 

or 

0*2  *  !  '*3a  \i  °l  CD2  V§  (Pt-Pa) 
and  similarly 

^c2  ’  I  l3c  Ap1  °1  c02  (Pt  ’V 

Thu*,  the  flow  from  the  flattened  area  L3j  D1  (or  l3c  )  is 
simply  the  normal  discharge  affected  by  Pt-P,  (or  Pt-Pc) 
reduced  by  a  factor  of  2/3. 

Now,  consider  the  condition  In  which  the  perforated  area  makes 
up  a  fraction  of  the  flattened  area.  The  diagram  in  Figure  50 
show*  the  possible  arrangement*  of  perforated  area  to  flattened 
area. 
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Conditions  1,4,10,11,12,  and  13  are  modeled  precisely. 

The  area  relationships  for  conditions  2, 3, 5, 6, 7,8  and  9  are 
accurately  represented  but  no  adjustment  has  been  made  to  the 
triangular  pressure  gradient  for  flow  in  the  flattened  portion 
of  the  trunk. 


Flow  Area  Adjustments  for  End  Elements. 


Ec2  '  (RN  *  7  L3c)  /R0 
Ed  ’  (*«  -  7  Lhcl>  <\ 

where  Efll,  Ea2,  Ec2  and  Ecl  are  the  area  adjustment  factors  for 
areas  al,  a2,  c2  and  cl  respectively. 
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8.  After  perforation  area  is  calculated  the  model  Indexes  to  the  next 
element  and  repeats  all  calculations  from  paragraph  1,  Section  4.2.2  In 
this  test.  Once  all  elements  have  been  evaluated  the  trunk  and  cushion 
volume  rates  are  calculated  through  Integration  using  a  first  order  lag. 

»c  "  t  (v;  •<c» 

\  ’  T  <V?-V 

where: 

v?  ■  calculated  volumes 
Vc,  Vt  *  volume  state  variables 
Vc,  Vt  »  volume  rate  of  change 

T  *  time  constant 

9.  The  total  effective  areas  for  air  flow  from  the  trunk  to  the  cushion, 
trunk  to  the  atmosphere  and  cushion  to  the  atmosphere  are 
calculated. 

Arel1ef  *  f  {pt  "  pa) 

CAca  *  cDgap  Agap 

2 

^ta  *  C01  Ahal  *  1  CD2  Aha2  *  CDa  Arel1ef 

CAtc  *  CD1  Ahcl  +  T  C02  Ahc2  *  CATV 

where  Arglief  *  trunk  to  atmosphere  relief  valve  area  and 

CA^y  *  effective  area  of  trunk/cushlorl  vent. 

The  relief  valve  aria  Is  a  function  of  the  pressure  difference  P  -P  and 

W  a 

Is  specified  by  the  user. 
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10.  Mass  flow  rates  are  calculated  from  subroutine  FNF10W  using  the  Chester 
Smith  compressible  flow  function  given  the  upstream  and  downstream 
pressure,  temperature  and  effective  area. 
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11.  Trunk  and  cushion  pressure  rates  are  evaluated  from  the  following 
functions  derived  from  the  conservation  of  energy  of  the  trunk  and 
cushion  volumes. 
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c  c 


■c  —  —  ”  'c 

Where  k  is  the  polytropic  constant  for  air. 
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SECTION  VI 


FOSTER-MILLER  TRUNK  ANO  CUSHION  COMPONENT 

6.1  General  Model  Description 

Foster-Miller  Assoicates  Inc.  have  developed  as  part  of  a  NASA  sponsored 
project  a  computer  program  to  analyse  and  predict  the  behavior  of  an  ACLS. 
The  program  structure  has  been  modified  and  made  a  standard  component  of  the 
EASY  ACLS  component  library.  The  analytical  models  upon  which  the 
Foster-Miller  program  was  based  are  described  in  References  6  and  7. 


SECTION  VII 


AIR  BAG  COMPONENT 


7.1  General  Model  Description 

7.1.1  Air  Bag  Geometry 


The  model  represents  the  geometry  and  pneumatic  characteristics  of  two 
parallel  inelastic  air  bags  (Figure  56).  The  air  bags  are  described  by  a 
finite  number  of  elements  which  need  not  be  uniform  in  cross  section,  width, 
or  location  on  the  vehicle  (Figure  57).  The  air  bag  elements  have 
unconstrained  lateral  movement  subject  to  ground  friction  forces  in  the  y-axis 
and  are  represented  with  a  membrane  model  as  sh  vn  in  Figure  58.  The  air  bag 
physical  dimensions  are  defined  for  the  right  air  bag  only.  The  pair  of  air 
bags  which  make  up  the  total  vehicle  supension  system  are  assumed  to  be 
symmetrical  about  the  roll  axis.  It  should  be  noted,  however,  that  the 
excitation  and  the  dynamic  response  of  the  air  bag  system  is  not  restricted  to 
symmetric  operation.  Each  air  bag  has  a  built  in  pressure  relief  valve  and  a 
provision  for  air  supply  as  shown  in  Figure  59.  The  element  locations  are 
specified  in  the  air  bag  coordinate  system  by  the  position  of  the  air  bag 
inboard  attach  point  "a"  (Figure  60).  The  air  bag  axis  is  referenced  to  the 
vehicle  body  axis  as  shown  in  Figure  56.  Thus  the  entire  air  bag  assembly  can 
be  repositioned  on  the  vehicle  by  changing  the  value  of  two  parameters,  BST 
and  WIT. 


7.1.2  Pneumatic  System 


A  schematic  of  the  air  bag  pneumatic  system  is  shown 
Compressible  flow  equations  are  used  to  describe  the 
remaining  air  supply  system  not  shown,  such  as  bleed 
are  modeled  with  other  EASY  components.  Air  outflow 
perforations  and  the  bag  relief  valve. 


in  Figure  61. 
state  of  the  fluid.  The 
air,  ejectors  and  fans 
occurs  through  the  bag 
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Figure  58  Membrane  Model  (Loaded  Shape) 
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Figure  59  Air  Bag  Supply  and  Relief  Valve 
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where  *AR  *  Total  Bag  to  Atmosphere  Flow  -  Right  Side 
■  Bag  to  Atmosphere  Flow  Thru  Relief  Valve 

\UB  •  Bag  to  Atmosphere  Flow  Thru  Lubrication  Holes 

“tr  ■  Air  Supply  to  Bag 


Figure  61  Air  Bag  Pneumatic  System  Schematic-Right-Side 


7.1.3  Component  Structure 


A  simplified  flow  chart  of  the  air  bag  component  (AB)  Is  shown  In  Figure  62 
which  describes  the  basic  function  of  the  model.  Input  variables  required  are 
velocities  and  displacements  from  the  vehicle  rigid  body  component  and 
flow-rate  and  temperature  from  the  air  supply  components.  During  the  inital 
pass  through  the  component  (l.e.  when  T«0)  parameter  data  tables  are 
calculated  from  a  set  of  Implicit  equations  which  represent  the  air  bag 
section  properties.  Subsequent  evaluation  of  the  element  section  properties 
can  then  be  made  economically  by  interpolating  these  tabular  data.  A  "DO" 
loop  Is  entered  within  the  code  which  evaluates  the  individual  element 
properties  and  calculates  their  accumulative  effect  as  an  air  bag  suspension 
system.  Note  that  the  terrain  model  Is  called  from  the  air  bag  model  and  is 
not  an  Independent  EASY  component. 

Vehicle  forces  and  moments  are  calculated  as  well  as  pressure  and  volume  rates 
of  change  and  flow  rates  which  are  outputted  to  their  corresponding  components. 

7.2  Discussion  of  Model  Idealization 

A  detailed  flow  chart  of  the  air  bag  component  is  shown  in  Figure  63  and 
corresponds  to  the  simplified  chart  of  the  model  shown  in  Figure  62. 

7.2.1  Air  Bag  Section  Properties 

Each  element  of  the  air  bag  is  represented  by  a  cylindrical  inelastic 
membrane.  The  set  of  equations  describing  its  shape  are  implicit  and  require 
an  iterative  solution.  Solving  these  equations  in  this  manner  for  each  call 
from  the  EASY  integrator  would  be  very  time  consuming  and  impractical.  The 
approach  used  for  this  model  is  to  initially  solve  the  set  of  equations  (i.e., 
at  T*0)  for  n  values  of  Ut  and  and  sto11*8  The  parameter  data  in  arrays.  The 
section  properties  are  then  evaluated  from  the  arrays  using  table  look-up 
routines  during  the  model  execution  (see  Figures  38  and  39). 
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gure  62  Air  Bag  Component 
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Input  TiriiPlpt: 


U.  *.  H 
t.  AIT 
R.  9.  I 
«■  I.  P 


tnm 
ioo  r 

Caappnppt 


*TR*  TT» 

VT7L 


Pro# 

Air  Supply 
Ti— oiwit 


1  RIGHT  AIR  RAC 
Ift  LEFT  AIR  RAG 


IS(t)  IS  THE 
ELEKKT  PARAMETER  SET 


OCTERWNE  APPROX. 
LOCATION  OP  POINT  *t* 


MATRIX  A  IS  THE  SOOT 
TO  EARTH  EULER  ANGLE 
TRANSFORMATION 


POP  ROOT  AM  EARTH 
AXES  *1  IS  COHN. 
TERRAIN  MODEL  *1 
IS  UR 


Figure  63  Air  Bag  Component  "AB"  Flowchart 
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LOJ 


•HAT*  PAAAXETtJB 
A «  ROTATED 
COINCIDENT  NOTH  THE 
AC  AIIS 


lo u  • 

■  Zo/St 

T.* 

St 

E(Z 

4* 

St 

nz 

4* 

St 

z 

f(Z 

S  * 

St 

f(Z 

rt  *  at4  <W 

rtx  •  “U  ft 
^  '-"ty  *t 


St  *  S.t 4  *> 

Fi  *  S  *  S, 

S  *  S  *  *fy 

rz*Fz  VS 

VV(VVVr^t 

z  V  VWV 


»T  ■  VT  ♦  *• 


CALL  lt»fl(  ) 


roETEWMES 
PEREDRAT10A  AREAS 

l  Si*  Si 


1  Viiitf  ’  '  <'r»-V 

*T1I 

•  Kl  •  V 

=*u 

•  Mj  AAj  ♦  |  CC 

aa2 

*  :a*r»H«f 

Citt 

(»>„.  V. 

tt»* 

cau;™>a») 

*TI 

"aa 

Sitta  . 

■*^r) 

f  REPEATED  FOR  I 

pi*  “a**  S*  j 


Output  VirtiPlu: 

rV  fV  Tz  T. 

6  OOF 

T  T  T  Cawpofltnt 

'*•  V  4 

Sr  Sa-  vtt  S* 

To  Alp 

V  Si-  Si-  Si 

"aa-  “al 

Supply 

Cocbo  n*f«U 

Figure  63  Air  Bag  Component  MAB"  Flowchart  (Concluded) 


The  section  parameter  data  is  generated  by  subroutine  ICB  as  shown  in  Figure 
64.  The  model  will  generate  and  store  up  to  six  sets  of  parameter  data  each 
representing  an  element  of  size  and  shape  specified  by  the  user. 

7. 2. 1.1  Element  Model 

The  free-body  diagrams  of  the  air  bag  element  membrane  and  geometry  for  loaded 
and  free  shape  are  shown  in  Figure  65.  The  following  assumptions  were 
observed  in  the  derivation  of  the  governing  equations. 

1.  The  air  bag  element  behaves  as  an  Inelastic  membrane. 

2.  Air  bag  membrane  tension  is  constant  along  the  free 
sections  of  Its  minor  periphery.  Contributions  from 
lateral  bag-ground  friction  are  Included. 

3.  The  air  bag  membrane  is  constrained  by  vehicle  hard 
structure. 

4.  Perforation  nozzle  reactions  are  negligible. 

5.  The  pressures  on  both  sides  of  membrane  element  In 
contact  with  the  ground  are  equal  to  Pt. 

6.  The  portion  of  membrane  In  contact  with  the  ground  Is  a 
straight  line. 

7.  The  membrane  is  tangent  at  the  edge  of  ground  contact. 

The  following  derivation  Is  for  the  free  air  bag  element  shape  shown  in  Figure 
65a. 


LQ  ■  r*j  +  r^  where  LQ  ■  membrane  peripheral  length 

ZQ  ■  rfl-cosa^-B 

Zg  ■  r(l-cos*2) 

A  ■  rsln^j  rsineg 
Y  ■  rslne* 

0  C 
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\  SUBROUTINE  ICB  INSET,  NPTS.  A.  B.  10.  GA.  GB.  DhU.  / 
V  AI0.AY0.ALl.AL3.AAS,  BTO,  BUT,  Bt3.  BAS)  / 


PURPOSE  -  SOLVE  FOR  AIR  BAG  ELEMENT  PAAAMETEPIC  MTA  FOR  FREE 

AND  LOADED  AIR  BAG  ELEMENT  SHAPES  AT  INITIAL  CONDITIONS 

METHOD  -SOLVE  MEMBRANE  GEOMETRT  AND  FORCE  BALANCE  EQUATIONS 
FOR  EVEN  INCREMENTS  OF  ZO  AND  MUT  AND  STORES 
PARAMETER  VALUES  IN  TABLE  LOCK  UP  ARRAYS. 

LIMITATIONS  -  CALLED  ONLY  BY  THE  AIR  BAG  MODEL  COMPONENT  AD 


—CALL  SEQUENCE— 

— •OUTF'UTS— 

ARRAYS  Y0.Ll.L3.AS  FOR  VARIOUS  10  AND  MUY 

•—INPUT  SU¬ 
NSET,  NPTS,  A,  B.LO.GA.GB, AND  DHU  (SEE  AIR  BAG  NOHEN.) 


WRITTEN  IT  -  J.R.KILNER 


APRIL  10.1370 


DO  10  I-1.N3ET 


P(3)-lll) 

P  01)  -LO  (I) 
P(7)-.017US-GB(I) 
P  (8)  -.  0174S-GA  (I) 
P 19)  -SIN  IP  (7)) 

P  (10) -COS  IP  (7)) 
Pin) -SIN  IP  18)) 

P  (12) -COS  IP  (8)1 
X  (l)  -0.25-LQ  (I) 

X  (2)  >2. 0 
X  (3D -2.0 
JQ-0 


C  SOLVE  FOR  FREE  SHAPE  PARAMETER  DATA 

l 


FSB.P. . 0001.  IP,  JQ.R.RMS.AJ.Z) 
sou3m 


IF 

(RMS.LC,  .OOOtl 
w  GO  tO  15  . 


TE  (IHRITE.  IN)  RMS.NSET  j 


Flyuro  64  Subroutln*  "ICB"  Flowchart 
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Figure  64  Subroutine  "IC8"  Flowchart  (Continued) 
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25 


crnmt 

gTO *X  12)  «SIN  (P  IIU3)-P  118)  U2) 

K^xiS 

ID  «X  (1) «  (.3 «P  (13)  *.  25*3  IN  12.  «P  (13) )  -SIN  (P  (13) ) ) 

"aSS  rl  S-P  1141  v 25-SIN  (2.-P  U41 )  -SIN  IP  m*l )  5 

a  m*  a .  -cos  ip  uaj )  >  *?  is  j  *p  iv  1 »?  a  si  «p  aoi 

«  *. 5*  a  (2) «  a . -CC3  IP  U4) ) )  *f  18)  )«PU8)-P(12) 
I)WS*nX*aXX*P  16)  «P  (2)  *,  5«P  (2)  «P  (3) 


(?  STUPE  POSITIVE 

— r 


a  swats') 


Figure  64  Subroutine  "ICB"  Flowchart  (Continued) 


|  WRITE  I1WRITE.  60201  1BL1  U.K,  I)  ,K«1 ,  NPTS) 


HS3 


CONTINUE | 

I 


WHITE  (IWRITE.  5034) 


X 


•H  W  su  J-l.NPTSl 

r _ —  J  .  - _ 

[  WHITE  (IWRITE.  60201  !flL3  U.K.  n  .K-l.  NPTSJ 
-  r- - 


"~tlL 


CONTINUE 

I 


WITE  tlWHITE.  60U2) 

I 


DO  55  J-l.NPTS  | 


|  WHITE  (IWRITE. 60201  (BL3  U.K.  I)  .K-l.  NPTSI 
- ~ - 


ss 


CONTINUE 

I 


1  WRITE  (IWRITE,  60U01 


f _ I _ , 

- - |  DO  60  >l.NPTS  | 

f _ _  i  - _ 

|  HMTE  (IWRITE.  60221  WHS  U.K.  11  ,K-1.  NPTSI 


- [60 


CONTINUE 

f _ =C~X_ 

~Z--1  r 

-H  go  si  J-l.NPTS 


X 


MITE  1IMITE,  60223  (IRQ  (J.K,  I)  .K*1.NPTS) 

_L 


1K. 


CONTINUE 

T 


-JlO  CONTINUE 

r 

I  RETURN 


Figure  64  Subroutine  "IC8"  Flowchart  (Continued) 
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Figure  64  Subroutine  "ICB"  Flowchart  (Continued) 
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gjiflflUTIMF  irLSBCX.M.K.A.PI 


PURPOSE  -  EVALUATE  ELEMENT  GECMETRT  AND  FORCE  BALANCE  EQUATIONS 

-  FOR  LOADED  SHAPE  AIR  BAG  ELEMENT 

METHOD  -  EQUATIONS  DESCRIBE  INELASTIC  MEMBRANE  SUBJECT 

-  TO  UNIFORM  INTERNAL  PRESSURE  AND  GROUND  REACTION 

LIMITATIONS  -  USED  ONLY  BY  ROUTINE  ONWT  IN  SUBROUTINE  ICB 

FOR  INPUTS  AND  OUTPUTS:  SEE  SUBROUTINE  ICFS 
WRITTEN  BY  -  J.fl.KILNEfl  APRIL  10,1378 


A(D«X(n-XC2)-PlS)«XlSl 
A  (2)  *P  (4)  -X  (3)  -X  (4)  -X  (5) 

A0)«P(2)-xm  «siN(p(i3n-xt2)«siNa»ii4n-xta 

■  *e  US)  «P  (10)  *P  (16)  «P  (12) 

A  (4)  -P  (3)  -X  (1) « (l . -COS  (P  (131)  1  «0C  (2) « (1  .-COS  (P I14J1 ) 

■  *PtlS)*»(9)-P(16)*«»(ll) 

A® -P  (6) -X  (2) -  (l -COS  (P  (14)) ) *P  (l©*4» (til 
RETURN 


Figure  64  Subroutine  "IC8"  Flowchart  (Concluded) 
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The  Free  Shape  equation  summary  and  problem  statement  are: 


Solve  iteratively  the  following  three  equations. 


Input  constants:  LQ,  A,  B 
Unknowns:  r, 

Determine:  ZQ  3  r(l-cos<*2) 

Y0  »  r(sin*2) 

4  ■  44 

The  following  derivation  is  for  the  loaded  air  bag  element  shape  shown  in 
Figures  65b  and  65c.  From  the  membrane  free-body  diagram 

ZT^sin 2(Pt-Pa)r2sin -j*  2 

or  Tt2  ‘  (Wr2 

Similarly  .  (P^P,)^ 

From  the  flattened  portion  of  the  membrane: 


where  F^  *  Element  ground  reaction 
F^y  ■  Lateral  friction  force 
0  ■  Element  width 
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And  is  the  instantaneous  friction  coefficient  which  is  a  function  of 
lateral  velocity,  i.e.. 


*t 


f(V,  ) 

fcT 


Note  that: 

% 

T  *  L3(IW 


Therefore  summing  forces  on  the  flattened  portion  of  the  membrane  yields 


or 

(pt-pa>rrMt4<pt-pa>-(pt-pa)'-2“0 


From  the  membrane  geometry  diagram: 


L0  «  L1+L2+L3 


where  L.  *  h.  +  r. 


‘I 


,  *  r.  ♦,  ,1-1.2 
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L1  -  P1(1T+Yi) 


h.  if  >  0 
0  if  h.  i  0 

k 

r 

TT+  Yi  if  h^  >  0 
Li/ri  lf  hi  -  0 


ZQ  *  r^fl-cos^-B-h.  sinY  1 
Zo  *  r2^1“COS4,2^~h2  sinY  2 

A  ■  Pj  si n 1-h1  cosY  1  +  sin*  2  "  h2  cosY  9  +  L3 
Thus,  the  Loaded  Shape  equation  sunmary  and  problem  statement  are 


’1  * 


’i  * 


rl 

-  rZ  4  ‘  0 

Lo 

■  4  *  4  *  l3 

A 

■  r1 sin*  t  +  r?sin*?  +  L,  - 

h^osY  ,  -  h?cosY7 

B 

■  r^l-cos*,)  -  r2(l-cos*2) 

-  hjSinYj  +  h2sinY2 

Z 

0 

-  r2(l-cos*2)  -  h2sinY2 

Input  Constants:  Ut.  LQ,  A.  B,  Zq.^.Yj 
Unknowns:  r 2,  ,  L2,  L3 

(Determine:  Yq  ■  r2s1n*2  -  h2cosY2 

L1 

1-3 

Aj  *  Y“2,  *j,  *2»  ^2) 

The  air  bag  element  cross  sectional  area  can  be  explicitly  determined  from  the 
membrane  equation  set  solution. 

a  .  1  4  _2.1.  2 

A1  7*lrl  +  7*2r2 

A2  •  |sfn2^1r12  +  j sin2*2r2Z 

2  2 
-  r^sln^  -  r2  sin*  2 

A3  "  7[rl^’COS0l)  +  Zo  +  B]hicosYi 
+  2[r2(l-cos*2)  +  ZQ]h2cosY2 

A4  "  ZoA  +  7  AB 


« 


A$  -  Ai  >  A2  +  A3  ♦  A4 


The  component  equations  and  logic  can  be  best  understood  following  the 
sequence  of  events  within  the  component  computer  code.  As  shown  in  component 
flow  chart  (Figures  62  and  63)  the  calculations  of  model  constants  and  air  bag 
element  parameters  are  carried  out  if  T=0  and  the  input  data  differs  from  the 
previous  case.  -Area,  volume  and  force  parameters  are  then  set  to  zero  for  the 
sunmation  procedure  in  the  following  loop  in  which  the  element  equations  are 
evaluated  for  each  element  in  the  air  bag.  For  each  pass  through  this  loop 
the  following  operations  are  performed: 


1.  A  test  is  made  for  right  or  left  air  bag  evaluations. 


E  ■ 


1.  right  side 

< 

-1.  left  side 


2.  The  approximate  location  of  point  "t“  is  determined  in  the  body 
axis.  (See  Figures  60  and  66).  First  the  free-shape  Zo  and  Yo  are 
determined  from  membrane  parameter  data  arrays. 
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^ofsu  *  AZ0(IS)  where  IS  *  Element  parameter  set 

Z  ,  a  s  7 
ofs  i  40fsu 

Yofs  s  Si  AY0(M.IS) 

where  subscript  “fs"  denotes  free-shape  and  subscript 
"u"  denotes  unsealed  variables  and  S.  is  a  user 
supplied  scaling  value. 

Then : 


XBt  3  XBai 

YSt  *  <Y8ai  +  YofS>E 
ZBt  3  ZBai  +  Zofs 

3.  The  location  of  point  "t"  is  then  determined  in  earth 
axis. 


Et 

-EtJ 


X 

[-ALT 


M. 


XBt 

YBt 

Z8t 


where  matrix  A  is  the  Body-to-Earth  Euler  angle 
transformation  (see  Figure  67). 


r1  ' 
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lzej 


r 
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V 

LZEJ 
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“cos'i' 

-sin^ 

0  “ 

V 

'x1' 

sin<|> 

_  0 

cos^ 

0 

0 

1  J 

Y1 

*  »  [ah*]  ■ 

Y  ’  * 

ZE, 

~COS0 

0 

sin0 

xb1 

0 

1 

0 

« 

Y 1 

‘  *  [A0] 

Y‘ 

-sinQ 

0 

COS0_ 

J'i 

L  J 

'  1 

0 

0  ■ 

fx 

*b 

w 

0 

cos* 

-sin* 

< 

yb 

>  3  [a*] 

yb 

_  0 

sin* 

cos*_ 

Zo 

l 

w 
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and 


M  ■  M  M  M 


4.  The  terrain  height  at  location  "t"  is  evaluated  from  subroutine 
TERRA.  Note  that  for  the  Body  and  Earth  axes  positive  Z  is  down 
and  for  the  terrain  model  positive  Z  is  up  (see  Figure  68). 

The  terrain  model  provides  an  option  for  defining  the  surface 
profile  with  a  sinusoidal  function  or  with  user  supplied  table  of 
profile  elevations.  These  methods  have  the  following  features: 

a.  Profile  Defined  In  Tabular  Form  (see  Figure  69) 

The  profile  data  must  be  specified  at  even  increments  along  the 
surface.  Linear  interpolation  is  performed  between  points. 

The  elevation  outside  of  the  defined  region  takes  the  value  of 
the  nearest  point  (see  flow  chart  in  Figure  47). 

b.  Discrete  (1-cosine)  or  Continuous  Sinusoidal  Profile 

Any  number  of  sequential  (i-cosine)  bumps  or  dips  may  be 
specified.  The  profile  begins  at  the  Earth  axes  origin. 

A  continuous  sinusoid  is  represented  by  a  large  number 
of  (1-cosine)  bumps. 

Thus 

Z£g  ■  TERRA  (X£t) 

5.  A  test  is  made  to  determine  If  the  air  bag  element  Is  In  contact 
with  the  ground  (see  Figure  68). 
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PROFILE  DEFINED  IN  TABULAR  FORM 


Figure  69  Terrain  Geometry 


lf  Zgap  i  0  Frw  Shape  •'loment  calculations  are  performed.  If  not, 
Loaded  Shape  equations  are  evaluated. 


a.  Loaded  Shape  Calculations 

The  velocity  of  point  "tH  In  Earth  axis  Is  determined.  The 
velocity  components  In  the  ground  plane  are  rotated  coincident 
with  the  body  axis  to  facilitate  the  ground-trunk  friction 
calculation, 

Z8t  *  Z8a1  +  Zo 

Then: 


fi  1 

r  ' 

f ' 

"ft 

$ 

( 

u 

P 

:Et 

'  -  C*«  CA»]  • 

V 

■  ♦  [8]  - 

q 

ZEt 

V 

w 

r 

»  4 

>. 

b 

VE  •Vi£2"  * 

t  6  t  fet 
Kto"  f<vEt> 


where  the  function  f  Is  the  mu-velocity  relationship  for  the  element  and  Is 
defined  as  shown  below. 


Next,  the  instantaneous  coefficients  of  bag-ground  friction  are  determined. 


V 


tx 


»*to 


k 


The  scaled  air  bag  section  properties  are  evaluated  which  are  a  function  of 
the  relative  distance  between  vehicle  and  ground  (Z  )  and  lateral  ground 
traction  (  jj  ). 


Zou  "  VSi 

Yo  *  Si  Wty^ 

L1  *  Si  WnVt/*) 

L3  “  Si2f4(Zou, 

As-Si  VZou,*VE> 


It  should  be  noted  that  for  the  purpose  of  calculating  element  membrane 
geometry  vrfien  in  ground  contact  it  is  assumed  that  the  surface  terrain  under 
the  element  is  parallel  to  the  horizontal  air  bag  axis.  The  resulting  error 
Introduced  is  minor  for  vehicle  roll  angles  less  than  10°. 


The  element  ground  contact  reaction  is  determined. 


A 

f 


t 

t 


DiL3 

<WA 


t 
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and  then  the  element  ground  friction  force  components  in  the  x  and  y  body  axis. 


Ffx  '  -  "tx  FT 

Ffy  ’  -  “ty  ft 

The  element  damping  force  is  evaluated  from  the  following  equation. 


Fd  *  °MP  ZEtL3Di 


This  relation  for  element  damping  is  based  on  the  following  assumptions. 


1)  Air  bag  element  damping  is  due  to  energy  dissipated  by  the 
deformation  of  the  visco-elastic  bag  material  and  the  air 
bag-ground  scrubbing  due  to  vertical  motion  of  the  vehicle.  It  is 
assumed  that  this  damping  force  is  proportional  to  the  relative 
vertical  velocity  between  the  vehicle  hard  structure  and  the 
ground.  Since  ground  rate  of  change  at  point  t  is  not  available 
within  the  model  the  damping  is  assumed  to  be  proportional  to  the 
absolute  vertical  velocity  at  point  t'  (see  Figure  66). 


2)  It  is  assumed  that  the  damping  force  is  proportional  the  footprint 
length  of  the  air  bag  element.  This  assumption  can  be  supported  by 
the  argument  that  at  the  instant  of  intitial  trunk-ground  contact 
the  damping  force  must  be  zero.  Whereas  at  maximum  bag  deflection 
the  bag  strain  rate  and  ground  scrubbing  must  also  be  maximum  for 
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constant  values  of  vehicle  vertical  velocity.  An  appropriate 
measure  of  these  factors  is  Lj. 

3)  It  is  self  evident  that  the  damping  force  must  be  proportional  to 
the  width  of  the  element,  .  The  damping  constant  of 
proportionally. 


is  then  a  function  of  the  retarding  force  per  unit 
velocity  per  unit  bag  contact  area  and  is  therefore 
independent  of  bag  size. 


The  ground  reaction,  friction  and  damping  forces  for  the  element  are  summed 
with  the  vehicle  forces  and  moments. 


F  -  F  +  F. 

X  X  _fx 

Fy  '  Fy  *  Ffy 
F*  *  Fz  -  Ft  -  Fd 


8t 


(YBai  *  Yo  +  7  L  )E3 


Tx‘Tx 


(Ft  +  Fd>  Y 


Bt  "  Pfy  ZBt 


Ty  "  Ty  +  ^Ft  +  Fd)  ^Bt  +  Ffx  ZBt 
Tz  “  Tz  +  Tfy  XBt  *  Ffx  YBt 


b.  Free  Shape  Calculations 


Only  the  element  cross-sectional  area  is  required 
for  the  bag  volume  calculation  if  the  bag  element 
Is  not  in  contact  with  the  ground. 

2 

AS  *  S1  Aas(1’1,IS) 
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6)  The  elemental  bag  volume  is  then  sunined  to  the  total 


7)  A  test  is  made  to  determine  if  the  air  bag  element  has  lubrication 
holes.  If  so,  the  orifice  areas  contributing  to  bag-to-atmosphere 
flow  is  calculated.  The  flow  chart  for  the  equations  and  logic 
used  in  these  calculations  are  shown  in  Figure  70.  The  assumptions 
and  equations  for  the  air  bag  out  flow  model  are  presented  in  the 
following  text. 


Definition  of  air  bag  perforations: 


A 


Figure  70  Flowchart  of  Perforation  Area  Distribution 
Equations  and  Logic 
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1  ■  ith  element 

*"hi  *  width  of  perforated  area 

Di  *  length  of  element 

'■pi  *  Peripheral  distance  from  attach  point  (b) 
to  edge  of  perforated  area. 

Api  ■  Orifice  are  per  unit  area  of  perforated 
surface  (or  porosity) 

Assumptions: 


a.  Air  bag  out-flow  through  perforations  not  in  contact  with 
ground  is  a  function  of  Pt,  Pa,  Cgj  and  the  corresponding 
portion  of  orifice  area. 


b.  Outflow  from  the  flattened  portion  of  bag  element  is  split 
based  on  the  location  of  a  stagnation  pressure  and  an  assumed 
triangular  gradient  and  is  a  function  of  Pt,  V  CD2  and  the 
corresponding  portion  of  orifice  area. 

c.  The  stagnation  pressure  in  the  gap  is  assumed  to  be  equal  to 
the  air  bag  pressure. 

d.  The  gap  pressure  gradients  on  either  side  of  the  stagnation 
point  are  assumed  to  be  equal  (see  Figure  71),  i.e.,oa  3 
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Figure  71  Pneumatic  Outflow  and  Pressure  Gradient 
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Therefore : 


VPa 

*F7 


3a  nr 
r— »  or 

L3b 


■3a  u3b 


*  7  *L3 


Now,  consider  the  condition  in  which  the  perforated  area  makes  up  a 
fraction  of  the  flattened  area.  The  diagram  in  Figure  72  shows 
the  possible  arrangements  of  perforated  area  to  flattened  area. 


Conditions  1,  3  and  6  are  modeled  precisely.  The  area 
relationships  for  conditions  2,  4  and  5  are  accurately 
represented  but  no  adjustment  has  been  made  to  the  triangular 
pressure  gradient  for  flow  in  the  flattened  portion  of  the  air  bag. 


With  the  assumption  that  the  pressure  gradients  in  the  air 
bag-ground  gap  is  triagular,  outflow  can  be  adequately  modeled  by 
reducing  the  normal  discharge  affected  Dy  (P„-Pj  by  a  factor  of 

U  a 

2/3. 


8)  After  perforation  area  is  calculated  the  model  indexes  to  the  next 
element  and  repeats  all  calculations  from  paragraph  1,  Section 
5i-2.2  in  this  text.  Once  all  elements  have  been  evaluated  the  air 
bag  volume  rates  are  calculated  through  integration  using  a  first 
order  lag. 


vt  ■  f  K  -  v 

where 

V*  ■  calculated  volume 

Vt  *  volume  state  variables 
■  volume  rate  of  change 
T  ■  time  constant 
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9)  The  total  effective  areas  for  flow  from  the  bag  to  the  atmosphere 
are  calculated.  The  equations  for  the  right  air  bag  are  presented 
below  which  are  identical  to  those  for  the  left  air  bag. 

Velief  *  f(ptr  '  Pa} 

CAta  *  CD1  Ahl  +  1  CD2  Ah2  +  C0a  Arelief 


where  A  ^  =»  bag  to  atmosphere  relief  valve  area. 

The  relief  valve  area  is  a  function  of  the  pressure 
difference  P  -P  and  is  specified  by  the  user. 

Li  a 

10)  Mass  flow  rate  is  calulated  from  subroutine  FNFLOW  using  the 
Chester  Smith  compressible  flow  function  given  the  upstream  and 
downstream  pressure,  temperature  and  effective  area. 


Y+l  1/2 


W  *  60  CA 

a  P 


r  7  ~7_'i’i 

ta  pt  TT  -  <T^>  J 


11)  The  bag  pressure  rate  is  evaluated  from  the  following 
function  derived  from  the  conservation  of  energy  of  the 
bag  volume. 


<"tr  -  »ar>  *  k 


Vtr 


where  k  is  the  polytropic  constant  for  air. 
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SECTION  VIII 


ARRESTING  SYSTEM  COMPONENT 


8.1  General  Component  Description 

Component  “AS"  is  a  dynamic  model  of  a  "Water  Twister"  manufactured  by  the  All 
American  Engineering  Co.,  type  of  arresting  system  as  shown  in  Figure  73.  The 
system  Is  a  simple  water  brake  that  converts  kinetic  energy  to  heat  through 
turbulence.  The  brake  consists  of  a  fluid  filled  steel  casing,  with  internal 
stator  vanes,  which  houses  a  vaned  centrifugal  rotor  (Figure  74).  The  rotor 
Is  mounted  on  a  shaft  which  extends  out  of  the  top  of  the  casing.  A  storage 
reel  for  a  nylon  tape  purchase  element  is  mounted  on  and  splined  to  the  top 
end  of  the  rotor  shaft.  The  tape  is  wrapped  on  this  storage  reel,  layer  on 
layer,  forming  a  spiral  wrap.  Pulling  the  tape  off  the  repl  causes  the  shaft 
and  vaned  rotor  to  revolve  within  the  fluid  filled  casing,  creating  turbulence. 

The  mathematical  model  can  be  divided  into  six  sections  for  the  convenience  of 
discussion.  They  are: 

(1)  Cable  and  tape  geometry 

(2)  Hook  and  cable  engagement  logic 

(3)  Tape  payout  and  strain  equations 

(4)  Hook  and  cable  impact  load 

(5)  Hook  and  vehicle  body  forces 

(6)  Energy  absorber  dynamics 

The  major  assumptions  made  for  the  derivation  of  the  arresting  system 
equations  are: 

(1)  Cable  strain  is  negligible  compared  to  the  tape  strain. 
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Figure  73  Water  Twister  Arresting  System 
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(2)  Cable  dynamics  are  not  significant  to  the  overall  landing  response 
of  the  vehicle  and  are  therefore  neglected  except  for  the  initial 
cable-hook  impact  resulting  from  the  change  in  momentum  due  to  cable 
mass  pickup. 

(3)  Tape  mass  is  neglected  except  for  its  contribution  to  the  drum 
inertia. 

(4)  The  hook  is  constrained  to  pivot  about  the  hook  attach  point  in  the 
X-Z  plane  of  the  vehicle  body  axis. 

(5)  Gravity  will  dictate  the  pre-contact  position  of  the  hook  relative 
to  the  vehicle  subject  to  the  constraint  stated  in  (4). 

(6)  Hook  slippage  along  the  cable  is  assumed  to  be  zero. 

Thus,  an  exact  mathematical  representation  of  the  arresting  engine  dynamics, 
the  cable/tape  geometry,  and  the  tape  steady  state  elastic  behavior  are 
included  in  the  model  with  a  simplified  representation  of  cable/tape  dynamics. 

8.2  Mathematic  Model 

The  development  of  the  component  equations  can  be  best  described  by  following 
the  sequence  of  operations  in  the  component  subroutine.  These  equations  are 
Illustrated  in  the  arresting  system  component  flow  chart  shown  in  Figure  75. 
Input  variables  required  by  the  component  include  the  vehicle  CG  position  and 
velocity  In  earth  axis  (X,  Y,  ALT,  XD,  and  YD)  and  the  vehicle  Euler  angles 
(*,«,  and a).  Component  output  variables  are  right  and  left  drum  angular 
displacements  and  velocities  («YR,  TL,  «R,  and**L)  and  vehicle  body  forces  and 
moments  (FX,  FY,  FZ,  TX,  TY,  and  TZ) . 

8*2.1  Cable  and  Tape  Geometry 

The  Inltal  cable/tape  geometry  after  hook-cable  contact  has  been  established 
1$  dictated  by  a  transverse  cable  wave  called  the  kink  wave  as  shown  In  Figure 
78  for  positions  1,  2  and  the  right  side  of  position  3.  Once  the  wave  has 
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reached  the  runway  edge  sheave  no  further  kink  wave  reflections  are  considered 
as  shown  in  positions  3  right  and  postion  4. 

The  kink  wave  angle  is  calculated  in  subroutine  KINK  and  is  called  once  during 
the  first  pass  through  component  AS  (i.e.,  at  T*0  and  if  the  input  constants 
for  KINK  differ  from  the  previous  case). 

Right  and  left  kink  wave  angles  are  assumed  to  be  equal  which  implies  that 
vehicle  sideslip  at  initial  hook-cable  contact  is  zero.  The  kink  wave  angle 
can  be  determined  from  the  following  implicit  equation  developed  in  Reference 

8. 


2  2 
«(!+«)-(£)  +  (VS(I^)  -e-e0) 

where  =  velocity  of  arresting  hook 


*  longitudinal  stress  wave  velocity 
1  *3C  in  the  cable 


Ec  *  cable  modulus  of  elasticity 
pc  *  cable  weight  density 
‘  *  cable  strain 
*o  3  initial  cable  strain 


The  equation  is  solved  implicitly  fore.  The  kink  wave  angle  is  then 


To  facilitate  the  computation  of  the  cable/tape  geometry  three  coordinate 
systems  are  defined  (see  Figure  77).  They  are 

(1)  Body  axis.  Located  at  the  vehicle  CG. 

(2)  Earth  axis.  Located  mid-span  between  the  runway  edge 
sheaves  at  the  runway  surface  level. 
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(3)  Hook-cable  axis.  Located  mid-span  between  the  runway  edge  sheaves 
at  Inltal  cable  height  above  the  runway  surface. 

The  hook-cable  plane  Is  defined  by  H,  L,  and  R.  Based  on  the  hook  and  cable 
kinematic  freedoms  P  must  be  in  the  hook-cable  plane  also. 


The  8ody-to-Earth  Euler  angle  transformation  coefficients  are  required  to 
determine  the  position  of  H  in  Earth  axis. 


CA]  -  CA*]CA0][A4] 


where 


[AO] 


cos* 

sin* 

L  o 


-sin* 

cos* 

0 


0 

0 

1 


[Ae] 


“  cose 

o 

-sine 


0 

1 

0 
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0 

cose 


[A  *3 


0  0  - 
cos*  -sin* 
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*  * 
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Y 

.  *- 

Y 

*  +  [A]  " 

0 

Z 

eh 

-ALT 

hh 

.  . 

k  a 

A  test  is  performed  at  this  point  to  determine  if  the  hook  has  made  contact 
with  the  cable  i.e.. 


H 

H 


<  0  then  no  cable  hookup 
>  0  then  there  is  a  potential  cable  hookup 


The  angle  between  the  hook-cable  plane  and  the  runway  surface  is  determined  as 
well  as  the  location  of  H  in  the  hook-cable  axis. 


H 


The  angle  must  be  determined  to  locate  P  in  the  hook-cable  axis  (see  Figures 
78  and  79). 

therefore  let 

a  *  X-  tanij)  see  Figure  79 

c  ■  Z^/cos  0 

Then 

fj-  ■  sin  0 


see  Figure  80 


Figure  79  Vehicle  and  Hook  Geometry  (View  1) 
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T~ - 


Therefore 


tanf  *  ^  *  sinQ  tan<fr 

Then  k  *  Z.  tane  +  Xr  /cosip 
L  lH 

b  *  d  -  k  tan£ 

s  K'  cos  C  _ b' 

D  *  D  cos(i|>-0  *  cos^i  +  sini^tanC 


1  p  tdn(i 

AB  a  3  +  b  *  X^tan\|)  +  C0Slj,+  sin,|,tanc  C  1 

-  tanC(Z.tan0+  Xp  /cosij/)] 

L  n 

g  » 

CH 

Finally 

sinfi  =  AB/g 

cos  5  3  X  /g 
lH 

The  location  of  P  in  the  hook-cable  axis  is  then 


see  Figure  81 


Tp  ricH 


8.2.2  Hook  and  Cable  Engagement  Logic 


A  test  is  performed  during  the  first  call  to  component  AS  in  which 

Xr  >  0  to  determine  if  the  hook  will  engage  the  cable  as  shown  in 
tH  “ 

Figure  82,  position  (1).  If  Xc  <0  then  engagement  is 

P 

assumed.  If  X  >0  then  the  hook  (point  p)  will  miss  the  cable 
P 

In  which  case  an  appropriate  diagnostic  statement  is  printed. 


202 


When  engagement  is  achieved  the  span  wise  intercept  constants  are  evaluated 
(see  Figure  78). 


h  ■  7  ys  '  *c 

\  •  7  »s  *  Tcp 

8.2.3  Tape  Payout  and  Strain  Equations 

The  cable/tape  lengths  from  L  to  P  and  R  to  P  irv  Figure  78  are 
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i-COS^y 

Al  '  TTrTJ^  »c.  *  j’s 


1 


'KW  -p 

The  tape  strain  for  the  right  and  left  sides  are 


*«■**»"  *  V 


*1  ■  <\  -  AIU)/(AIU  *  V 

where  AR(J  *  unstrained  tape  payout,  right  side 
Aq|  »  unstrained  tape  payout,  left  side 

Tape/cable  tensile  loads  are  then 


TR  ■  1  “  t  (CR> 

TL  '  1  “  t  <V 


where  t  *  tape  thickness 
W  «  tape  width 

Et(e)  »  tape  stress  as  a  function  of  tape  strain 
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8.2.4  Hook-Cable  Impact  Load 

When  the  vehicle  arresting  hook  strikes  the  cable  a  kink  wave  is  generated  as 
shown  in  Figure  76.  The  cable  impact  load  resulting  from  the  cable  dynamic 
response  is  approximated  by  the  change  in  momentum  due  to  the  cable  mass 
pickup  during  the  period  between  first  contact  and  the  first  reflection  of  the 


V  «  VxD2  +  YD2 
s 

and  that  n  *  YD/XD 
for  the  left  side 


where  FL  *  impact  load,  left  side 

■  rate  of  cable  mass  pickup  by  the  hook 

Then 

cfcn  »  Ac  pc  da 

dt  g  dt 

where  Ac  *  cable  cross-section  area 

pc  »  cable  weight  density 

Note  that 
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sin(n/Z-$KW-n)  sin^KW 
cos(^  +n) 

or  a  *  — -t - b 

S1n*KW 

da  *  cos(^u  +  n )  db 

and  -rr  - = — - V  ,  where  V  *  -rr 

dt  sin  ♦  KW  s  s  dt 

Therefore 

Acpc  cos(*in,+  n)  ? 

F,  * - - — - V  c  if  X_  >  0  and  Avu>  6  , 

L  g  Sin  ^KW  s  Eh  rKVI  r  L 

Acpccos(4  -n)  2 

F„  - - : - V  *  if  X_  >  0  and 6  u>  0  0 

R  g  sin$KW  s  rKW  *  R 

8.2.5  Hook  and  Vehicle  Body  Forces 

The  hook  angle  relative  to  the  body  axis0H  is  required  to  resolve  the  hook 
forces  into  body  forces  applied  at  the  vehicle  CG  (see  Figures  83  and  84). 

Let  Z«  equal  an  arbitrary  constant  length  parallel  to  the 
dG 

Z  body  axis  from  H  to  G.  Then  the  location  of  G  in  the  Earth  axis 
Is  _  _ 


♦  CA]  J  0 


1  eg  z  eh 


The  location  of  I  in  the  Earth  axis  is 


Y  ,  *.  Y  c-(a+b) 

Z.  E  I  “H  C 


Y2  •  (XE  -XE  )2  ♦  (Ye  -Ye  )2  ♦  (ZE  -ZE  )2 

eg  ei  eg  ei  eg  E1 
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Figure  84  Hook-Body  Axis  Angle 
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8.2.6  Energy  Absorber  Dynamics 

The  equations  for  the  energy  absorber  are  derived  for  the  unit  on  the  left 
side  of  the  runway  and  are  identical  for  the  right  side  (see  Figure  85).  The 
equations  of  motion  are 

<!  *  W“l  ■  -  °MP  "L  I"lI  *  *1.  TL 
tl  *"l 

where 


j 


\ 


^DR  *  Inerti a  of  the  rotating  mass  except  the  tape 
DMp  »  Oamping  coefficient 
R^  *  Instantaneous  radius  of  stored  tape 
T^  »  Tape  tension 

»  Angular  velocity  of  the  drum 
*  Angular  displacement  of  the  drum 

The  contribution  of  the  stored  tape  to  the  total  inertia  is 


JHEL  4_r 
7  g  lKt  KI  > 
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P 


'i*  .yrmvm 


rrr 


where 


W  *  Tape  width 
Oy  *  Tape  weight  density 
Rj  *  Radius  of  Tape  hub 

The  tape  radius  is 


and  the  unstrained  tape  payout  is 


LU 


"oV1’ 


T  L 
47TR 


where 


t  ■  tape  thickness 

R„  *  Drum  outside  radius 
o 
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